
DELHI UNIVERSITY LIBRARY SYSTEM 


GS 


Gi 


Cl No 

Ac No cfO <^S ^ Date of release of loan 

This book should be returned on or before the date last stamped 
below An overdue charge of 20 paise will be charged for each day the 
book IS kept overtime 




CARTESIAN TENSORS 



LONDON 

Cambndge University Press 
PETTEB LANE 

NEW YOBK^TOBONTO 
BOMBAY • CALCUTTA • MADRAS 
Macmiiian 

TOKYO 

Maruzen Company Ltd 
AU nghta reserved 



CARTESIAN TENSORS. 


BY 

HAROLD JEFFREYS 
M A, D Sc. F R S 


CAMBRIDGE 


AT THE UNIVERSITY PRESS 
1931 



PBINTED IN QBEAT BRITAIN 



PREFACE 


It IS widely felt that when the equations of mathematical 
physics are written out in full Cartesian form the structural 
simphcity of the formulae is often hidden by the mechamcal 
labour of writing out every term explicitly Attempts have 
been made to reduce this labour by one form or another of 
vector algebra, but it has always seemed to me that this 
method both mtroduces new difficulties and is insufficiently 
general Thus the product of two vectors, in vector lan- 
guage, means one of two thmgs, either the scalar or the 
vector product, and it is not physically obvious why just 
these functions of the vectors should anse and no others 

The use of tensor notation, with the summation con- 
vention, carries out as great a simphfication of the wnting 
as does vector notation The notation has actually at- 
tracted attention owing to its applications in the theory of 
relativity, but for ordinary purposes two great abbrevia- 
tions may be made We use rectangular Cartesian axes, 
the result is that the distinction between covanant and 
contravanant vectors disappears, and with it the terms 
arising from curvature of the surfaces of reference The 
formidable character of most of the formulae of the theory 
of relativity is absent from the formulae of tensors referred 
to Cartesian axes The tensor method is a necessity for re- 
lativity, for applications in dynamics, electricity, elas- 
ticity, and hydrodynamics it is a great convemence. 

It IS found that the scalar and vector products are not 
the only functions of two vectors that arise, though the 
theory provides reasons why they are important m many 
apphcations There is also a symmetrical product, which 
ordinary vector notation is completely unable to express. 
In tensor notation it arises naturally as a symmetrical 
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tensor of the second order. The system of moments and 
products of inertia of a ngid body constitutes such a 
tensor, so do the stress components and the strain com- 
ponents in an elastic sobd. 

The present method, bke vector notation, is of use prin- 
cipally in proving general theorems In concrete apphca- 
tions there is usually some asymmetry about the coor- 
dinates that makes it necessary to abandon the tensor 
form at some stage in the work It has been said that 
vector equations are hke a pocket map, and it has been re- 
phed that a pocket map has to be taken out of the pocket 
and unfolded before it is of any use. The same apphes to the 
tensor method, and for the same reason, but it has the 
great advantage that it is not a new notation, but a concise 
way of writing the ordinary notation, so that the unfolding 
can be earned out more convemently when occasion arises 

What 18 usually called Statics is treated in Chapter v, 
after Dynamics I consider this to be the proper order, be- 
cause Statics IS a special case of Dynamics, and many of its 
formulae have physical significance for reasons explained 
in Dynamics The customary reversal of the order is due, 
I beheve, to the fact that an introduction to mechames has 
to be given at schools before the students have received 
any traimng in calculus, but this need not influence 
students working for a umversity examination 

It should perhaps be stated that the object of this work 
IS to illustrate the use of tensor methods , it does not claim 
to give a complete theory of all the subjects touched, re- 
ference for which must be made to the standard text-books 

I must express my gratitude to Mr M H A Newman, 
Miss L M Swam, Dr S Goldstein, and Dr Bertha Swirles 
for assistance at various stages in the work, and to the 
staff of the University Press for their care m the prmtmg. 

HAROLD JEFFREYS 

September 1931 
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CHAPTER I 


CARTESIAN TENSORS 

If we have two sets of rectangular axes (Ox, Oy, Oz), 
(Ox' , Oy', Oz') at the same origin, the coordmates of a 
point P with respect to the second set are given m terms 
of the coordinates with respect to the first set by the 
equations 

x' = l^x + rriiy + tiiZ' 

y' = l^x + m^y + n^z ■ (1) 

z' = igx + m^y + WgZ 

The quantities (li,mi,n^, ..,n^) are the cosines of the 
angles between the vanous axes, thus is the cosme of the 
angle between the axes Ox' and Ox, is the cosme of the 
angle between Oy' and Oz, and so on It follows that the 
coordinates (x, y, z) can be expressed in terms of (x', y', z') 
by the relations 

X = l^x' + l^y' + I3Z' 1 
y = rriix' + m^y' + m^z' ^ ( 2 ) 

z = n-^x' + n^y' + n^z' ] 

We can shorten the wnting of (1) and (2) considerably 
by a change of notation Instead of (x, y, z) let us wnte 
(Xj, X 2 , X 3 ), and instead of (x', y', z') wnte (x^, x^', x^) We 
can now say that the coordmates with respect to the first 
set of axes are x,, where i may be 1, 2, or 3, and those with 
respect to the second set are %! , where j may be 1, 2, or 3 
Then m ( 1 ) each coordinate Xj is expressed as the sum of 
three terms depending on the three x, . Each x, is associated 
with the cosme of the angle between the direction of that 


J C T 


I 



j2 CARTESIAN TENSORS 

Xf increasing and that of x/ increasing Let us denote this 
cosine by . Then we have, for all values of J, 

Xj ~ ^ 31^3 

= S a^jX, ■ (3) 

i-l. 2,3 

Conversely (2) can be written 

x^= L a,jx/ (4), 

J = l. 2,3 

the a,j having the same value as in (3), for the same values 
of t and j, because it is m both cases the cosine of the angle 
between the directions of x, and x/ increasing. 

In mathematical physics we often have to deal with sets 
of three quantities in relation to a set of axes, of the general 
form u, (that is, u^), and such that in relation to a 

different set of axes the corresponding quantities are 
(iq', ttj', Mj'), which satisfy the relations 



It/ = 2 a,jU, 

(5) 


z=l,2, 3 


and 

It, = 2 a„u/ 

(6) 


1 = 1. 2.3 


Such sets of three quantities are called tensors of the first 
or^r, or vectors The mdividual %, it, may be called the 
components of the tensors 

Clearly if we multiply all of the it, and uf by the same 
quantity m we get 

mu,' = S [mu,) (7), 

1=1. 2,3 

so that mu, is another tensor of the first order 

Agam, if we have two tensors of the first order, it, and 
Vi , we shall have 

«/ + = S a,i (it, + Vi) (8), 

z-l, 2, 3 

so^at it( + is a tensor of the first order. 

'^We notice that each of the equations (3) to (8) is really 
a set of three equations, where the suffix i otj appears on 
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the left it is to be given in turn aU the values 1, 2, 3, and tlfe 
resulting equation is asserted in each case. In each such 
equation the nght side is the sum of three terms, obtained 
by giving j or ^ the values 1, 2, 3 m turn and adding Wher- 
ever such a summation occurs a suf&s: is repeated in the 
expression for the general term,, where there is a summa- 
tion for all values ofj the general term, such as atjU/, con- 
tains j twice We make it a regular convention that, unless 
the contrary is stated, (whenever a suffix is repeated it is 
to be given all possible valu^ and that the terms are to be 
added for all^ Thus we write (5) as simply 

< = ( 9 ). 

the summation sign being automatically understood by 
our convention Then (9) really means three equations, 
with three terms on the right of each, but we can by means 
of our conventions express all of the twelve terms com- 
pactly by the single equation (9) 

There are single quantities, such as mass and distance, 
that are the same for all sets of axes These are called 
tensors of zero order, or scalars. 

Consider now two tensors of the first order, w, and v,^ 
(When we write “ a tensor u,” we mean of course a tensor 
of the first order whose components are This is 

another piece of shorthand.) Suppose each component of 
the one multiphed by each component of the other; then 
we obtain a set of mne quantities expressed by u, , where 
each of i and k is independently given aU the values 1, 2, 3 
The components of u , , f * with respect to the other set of 
axes are uf, Vi' say , and 

« = (OfyM.) (OtiWjfc) I-- '''" ^ , 

= (19) 

The suffixes i and k are repeated on the right. Thus (10) 
represents mne equations, each with mne terms on the 
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rtght. Each term on the right is the product of two factors, 
one ^ the form , dependmg only on the orientation of 
the axes, and the other of the form u^v^, representing the 
products of the components referred’ to the original axes. 
In this way the various Uj'vi' can be obtained m terms of 
the ongmal u, . But products of two vectors are far from 
being the only quantities satisfymg this rule. In general a 
set of mne quantities referred to a set of axes, and 
transformed to another set by the rule 

fll), 

IS 6auld a tensor of the second order. 

We may go on similarly to construct and define tensors 
of the third, fourth, and higher orders Thus a set of quan- 
tities that transforms like x^x^x^Xj, is called a tensor of 
order n, where n is the number of factors in this product 

When we say that a certain set of quantities is a tensor 
of any order n, we mean that we have ways of specif3ang 
its components with respect to any set of axes, and that the 
components with regard to any two different sets of axes 
are related according to the rule appropriate to tensors of 
that order, and in particular to the products of the coor- 
dinates with n factors For instance, if we say that m, is a 
tensor of order 1, we are not simply defimng uf as meamng 
a,>u,. We are supposing both that u/ has a meamng, such 
as a displacement or a velocity, with reference to the axes 
of xf , and that the value of each component is equal to 
Thus the statement that any set of quantities is a 
tensor is not a mere convention, but a statement capable 
of test and therefore needing proof In (7) and (8), for 
example, our data are that u, and m/ are the components 
of a vector with regard to two different sets of axes We 
prove that the sets of quantities obtained by multiplying 
both by the same quantity are related accordmg to the 
vector rule, and therefore the products are vectors. 
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If we interchange^' and I in (11), we get • 

Wi/ — aiiOrfcjWi^ ( 12 ). 

But on the right i and h are “dummy suffixes”; that is, 
they are to be given all possible values and the results 
added. It is unimportant which of them we call % and which 
k ; we may therefore interchange them and get 


M5[,' = — ayajktW>)|., (13). 


Thus transforms accordmg to the same rule as and 
therefore is another tensor of the second order The im- 
portance of this 18 that if we know the set of quantities 
arranged 


/M’U 


W’ls' 

( “>21 


“’23 

Wai 

M’sa 

W33, 

the second order. 

/M'll 

Wa 


1 ^12 


W'sa 





(14) 


(15) 


IS another tensor of the second order Therefore the sets 
{w,k -t- w,c,) and {iv,^ — «;*<) are tensors of the second order. 
The first of these has the property that it is unaltered by 
mterchangmg i and k, and is therefore called a symm etricgJ, ^ 
tensor The second has aU its components reversed m sign 
when i and k are interchanged, and is called an avMs ym - 
metncal tensor Clearly in an antisymmetncal tensor the 
“leading diagonal” components, le those with i and k 
equal, are all zero Also, since 

Wit = i {Wii -f- Wfc,) 4- i (wjfc- wt,) (16), 


we can consider any tensor of the second order as the sum 
of 83nnmetncal and antisymmetncal parts 
The gradient of a scalar is a vector. For if is a scalar, 
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itfe gradient is dU/dxt or dUjdXf' according to the set of 
axes/ But 


dx, w 

dx/ dx/ dx, dx. 


fso that the gradients transform according to the vector 
^rule. Similarly the gradient of a vector is a tensor of 
^ order 2 For if u, and %' are the components of a vector 
with respect to two sets of axes, 


du/ 

dx/ 


dx^ du/ d , , 

d^/d^~ “*' 0 ^ 


_ du, 

~ dx. 


( 18 ), 


so that the rule of transformation is as in (11) 

I Since X, IS a vector, it follows that dxjdx,. is a tensor of 
the second order But dx,ldx, is uiuty li i = k and zero if 
t ^ k Hence the set of quantities 8,,, such that 



constitutes a tensor of the second order We can prove this 
directly, for if we apply (11), 8,/ m the new system of 
coordinates should be given by 

^3l — ( 1 ®) 

The suffix k has to take all values 1, 2, 3 But if fc # i, 

8,, 18 0, and the corresponding term is zero 

8,;^ = 1, and the result of the summation with regard to k is 

)) 5 ,(' = ( 20 ) 

But the o„ are the direction cosines of the axis of x/ with 
regard to the x,, and the a,i are those of x/ with regard to 
X, Hence is the cosine of the angle between x/ and 
x/, and is equal to 1 if the axes are identical and to 0 if they 
are perpendicular It follows that the result of the trans- 
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formation is that 8^’ = 1 if j = Z, and 8^' = 0 ii j ¥= I- Jt 
follows that the set of quantities 



IS transformed into itself by the rule (11) and therefore is 
a tensor of the second order 

If Ui is a vector and we form the product have 

a tensor of order 3. But now put m = k and add for all 
values of k Smce 8.^ = 0 except for k = i, the only term 
different from zero is that for k = i, and this is . Hence 

= u, ( 22 ) 

This operation therefore replaces the suffix k by i The 
tensor S,* can therefore be called the subshtuHon tensor 
^ In the Tensor u;,*, let us put k= %, an^ in 'accordance 
with our convention add for all values of k. Then the 
corresponding quantity w,/ is got by putting I = j and 
summing, but 

w„' = 

= M’.. (23) 

Thus w„ transforms into itself and therefore is a scalar^J 

This operation of putting two suffixes in a tensor equal 
and adding accordingly is known as contraction In general 
it gives a new tensor, whose order is less 2 than that jji 
the origmal tensor If for instance we contract the tensor 
M.Wj., we obtain 

u,v^ = u^Vi+ + UgVj (24), 

which is the scalar product of tt, and 

Similarly the same tensor yields the symmetrical and 
antisymmetrioal tensors («<«*, + u^v,) and («<«* — 

We may call these the symmetrical and antisymmetrioal 
products of U( and v* . 
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The t^psor Sujt/dx, gives similaxly, on contraction, a scalar 

p»), 


dx, dx^ ^ 0Xj ^ dxg 
which is known as the divergence of while it gives also 
symmetrical and antisymmetrical tensors 

0% du, , 0Mj. duf 
dx, 0a;* dx, 0a;* ' 

The former has important apphcations, especially in the 
theory of elasticity and hydrod 5 mamics , the latter is 
known as the curl or rotation of u. The vamshing of the 
curl IS the condition that u, may be the gradient of a scalar 
All the above considerations can be extended to any 
number of dimensions In n dim ensions a tensor of order r 
has components A tensor of order 2, m particular, has 
components If it is antisymmetrical, the n diagonal 
components are zero, and the others are equal and opposite 
m pairs Hence an anti 83 nnmetncal tensor of order 2 has 
(» — 1) independent components If w = 1, 2, 3, 4, , in 
turn, this number is 0, 1, 3, 6,... It happens that m three 
1 dimensions the number of numencally independent com- 
ponents of an antisymmetrical tensor of the second order is 
equal to the number of components of a vector Actually it 
can be proved that with any vector we can associate an anti- 
symmetncal tensor of the second order, and conversely 
This IS not true in any number of dimensions other than 3. 

Smce the a„ are the direction cosines with respect to the a;< 
of three perpendicular hnes, they are connected by six re- 

au* + + 0*1=' = I 

«12* + ^ ^ (26). 

« 13 * + “ 28 * + 0 * 8 * = 1 

0i 20*3 -f- ^22^23 “f~ Cl32®83 ~ 0 

®isOii -l- ®2s02i "f" OasOsj =0 ■ (27). 

On Oi2 ~ f ~ ^ 21 022 O31 032 ~ 6 ^ 


0a* 0M, 

dx, 0a;* 
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We notice that the second and third of (27) both*contain 
(«U 5 <hi)- We may therefore solve them for the ratios 

of these quantities. Thus 


^21 


^ _ jj. 

^ 38^2 — ^ 8®82 ® 13^22 — ® 23®12 


say Substituting in the first of (26) we get 


(28), 


1 = 

^ {®U {®28 

®32 ~ 

Clg^C 

= 

- * 

®11 

®21 

®31 



®12 

®22 

®32 

Also 


«13 

®28 

®28 

{(^23^^22 “ (®33®12 


+ (<^ 18^22 ~ ® 23 ® 12 )} 

(29) 


(30) 

But we have a general identity 

(a^ + 6* + c*) (o'* + + c'*) - (aa' + bb' + ccy 

= (be' - c6')" + (ca' - acy + (ab' - ba')^ (31) 

Hence 


** [(®13* + ®23* + *88*) (“m** + »28** + «32®) 

(®ia®13 “f" ®22®23 ®32®ss)*] ~ 1 (^2) 

But on account of the second and third of (26) and the first 
of (27) the expression in brackets is umty, and therefore 

* = ± 1 (33) 

For any given transformation the determinant in (29) is 
therefore equal to ± 1. Evidently its sign is reversed if we 
interchange any two of the suffixes for this mterchanges 
two rows of the determmant^so that the sign is a matter of 
the numbering of the axes ,If we start with a rigid frame 
attached to the axes x,, and rotate it contmuously till it 
18 attached to the axes x/, all the o<y vary contmuously and 
therefore the determinant cannot change from + 1 to — 1 
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or from — 1 to + 1 If then goes to Xj", x^ to x^, and ajg 
to x^, the determinant is initiaDy 


and therefore 


1 0 0 
0 10 
0 0 1 


1 


(34), 


I; = - 1 (35), 

and the determinant formed by the a,, is always + 1 
If we have a frame of axes (XiX^x^) we can turn it by a 
continuous movement so as to bring Xj along the old Xg, 
X 2 along the old x ^ , and Xg along the old Xj In this case we 
have 


and 

II 

II 

II 

(36), 

®ii = 0, 

®21 — 1, ®31 — ^ > ®12 — ®22 — ®32 — ^ > 

®13 ~ 1) ®23 ~ ®33 ~ ^ 

((37) , 


/ ■ 


The determinant of the a„ is therefore 

I 0 1 0 I = 1 (38), 


0 0 1 


1 0 0 


as before Any rotation of the axes that does not alter the 
cychc interchange of suffixes 1 2 3 1 2 3 1 . . therefore 
leaves the determmant equal to imity, and therefore so 
long as we always use right-handed or always left-handed 
axes the determinant of the is - 1 - 1. 

With this restriction 


— ®22®33 ®32®23> ®21 — ®13®32 ®S3®12 

^31 “ ^ 23^2 — 



and therefore every direction cosine is equal to its first 
minor in the determinant 
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These relations are of course identical with those ex- 
pressed in the usual notation of sohd geometry by 

h = — n^li, % = hni^ — 

Now suppose that is a vector, and consider the set of 
quantities 

«>.*=/ 0 M 3 -u^ (41), 

I - M3 0 Ml j 

\ Mj — Mj 0 / 

Apply (11) to this, taking i to be the number of the row and 
k that of the column We see that enters as and as 
— MJ 32 . Its coefficient in m>,,' is therefore 02,031 — 03 , 021 , 
and in all 


'^ll — (®2J®3! — ®2l) ^1 + (®31®1! ~ ^2 

+ (01,021- «2J«1!)«8 (42). 

This is obviously zero if j i and if the other axis 

perpendicular to x,' and x/ is x„', and J In jin is a cychc 
order, the quantities in brackets are equal to (Oi„ , Ojn , 03 ™) 
This 18 true if I immediately succeeds j in the order If I 
precedes j by one place the signs are reversed Hence if 
j = 1 and Z = 2, or J = 2 and Z = 3, or if j = 3 and Z = 1, 

jj w,i = o,„M, = M„' -i (43), 

and in the alternative case 

Wfi' = - M„' (44) 

Thus w,i' = / 0 Mj' — Mj'X (45), 

( - M 3 ' 0 uA 

\ U/ - Ml' 0 / 


and is of the same form as (41). Thus with any vector we tj 
can associate an antisymmetncal tensor of the second order. , 
i Conversely with any antisymmetncal tensor of the second 
order we can associate a vector 
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^ We can proceed alternatively by considering the set of 
quantities eami defined by the condition that if any two 
of t, k, m are equal the corresponding component is 0; if 
i, k, m are all unequal and m cychc order, the component 
IS + 1 ; if the order is not cyclic, the component is — 
Let us see whether this is a tensor of the third order. If so, 
we should have 


— 02/^21^31 + ®2l®Sl®ln + 

— (hjO^idan ~ ~ ® 1 #® 8 !® 2 « (^®) 

Now if, for instance, j = I, the nght side is clearly zero and 
^jin' = 0 If j, I, n are all unequal, the expression is 


®iy ®2y ®8y 

Oil Ugj Ojj (47), 

®ln ®2n ®3n 

which is equal to 1 if jin are in cychc order and to — 1 if 
not. Hence the set of quantities is transformed mto 
* itself by the rule for transforming tensors of order 3, and 
therefore constitutes a tensor of order 3 This is called the 



alternating tensor. 

Now consider the product where is a vector 

This is a tensor of the fourth order. If we contract it by 
puttmg p = m and summmg we get a second order tensor 
W’i* = If » = 1 and k = 2, the only value of m that 

makes difierent from zero is 3, and then = + 1 Hence 

«'i2 = (48) 

If ^ = 2 and k = 1, m is 3, but 213 is the reverse of cyclic 
order and e„a = — I Hence 


— “^3 


(49). 


Similarly we find that the elements of are 
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ao that the antisymmetrical tensor aaaociated with a vector 
jcan actually be obtained from it by multipl3ring by 
|md contracting. - ' ' , - ,, 

y Again, suppose that we are given a tensor of the second 
order and that we form a vector Ut by multiplying by 
and oontractmg twice. We have, if wi = 3, 


■*^3 — ^ika'^tk — ^ 128*^13 + ^218 ^21 — *^12 — 

Thus (61). 

If Wfjc 18 symmetrical, this evidently gives zero. If it is 
antisymmetrical the components of are numerically 
twice those of 10 ^^. 

On account of the intimate relation between the vector 
and the antisymmetrical tensor we shall habitually denote 
the tensor w,^. of (41) by m,*, so that 

UiJ = ligj = 1133 = 0 , lijj = Ms , 'U33 = tij , llju = , 

^21 = W 3 , ~ — , "Mja = ^2 

It will always be seen at once whether the vector or the 
tensor is mtended, smce the former has one and the latter 
two suffixes. 

If we have any three vectors u^,v,,w^, and consider the 
scalar we see that 



1’ = M1V2W3 + U^V^Wi + 



— 


- UjV^Wi - U^VsWi 


II 2 

«3 1 


(53) 



^3 



M’l 

W 2 

«’$ 1 




so that we have a concise way of wntmg the deter minan t 
formed by the components of three vectors If any two of 
the vectors are parallel this scalarjvaniahes. 

In associatmg a vector with an antisymmetrical tensor 
of order 2 a sign convention clearly arises. We make the 
positive signs m (41) he one place to the right of the leading 
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iagonalyTLf then we have two vectors Ut and v^, their anti- 
symmetrical product IS M.v*. — u^Vi, and m the associated 
vector we give the positive sign to when k follows i in 
the cyclic order. Hence the components of this vector are 
taken to be 


(■UjjWa - «3t;g, - UyV^, UiVi - ttgVi) 


This vector is perpendicular to both the original vectors , for 

Ml (UiVs — V^V^) + (M3M1 — Mi^s) + W3 {u^V^ — = 0, 

Vi (MjjWj - v^Vi) + 1>2 {u^Vi — MiVg) + Vg (Mi1;jj - u^Vi) = 0 


We call it the vector prodvict of «, and Vu> and can save 
writing by denoting it by [«, «]„ 


3mj, 


we 


Similarly xvith the antisymmetncal tensor ^ 

CX^ OXfc 

associate a vector so as to leave the sign unaltered when 
k follows % m the cyclic order Thus the components are 


/dv^ _ Smj dv^ dug dug Smj 
\ 3 x 2 dza ’ 8X3 8xi ’ 0Xi 8X3 


(54) 


This is often called “curl u ” 

The Tensor eaafmvt- Since this tensor is the product of 
two third order tensors, once summed, it is a tensor of 
the fourth order, i, k, m, p being arbitrarily assignable 
Evidently li i^k ot m = p, the corresponding component 
IS zero 

If i = m, the contribution from any value of a is zero 
unless also k — p, and then 


and the component is + 1. 

If i = p, then no value of s gives a contribution unless 
k ~ m. Then one of e,*, and is + 1 and the other — 1, 
and the component is — 1 Hence the components of the 
tensor are as follows 



CARTESIAN TENSORS 

li i = m,k = p, the component is + 1, unless i also = 
i p, k = m, the component is — 1, unless » also = 
i = k OT m = p, the component is 0 
These results apply also to the tensor 



k, 


and therefore 


^tka^mpa “ ^im^kp ^tp^km 


(55) 


EXAMPLES 


1 If M,, u/, Uj^’ are the coiupoaenta of a vector with regard to 
three sets of axes, prove that the values of are the same as 
would be obtamed by transforming first from % to u/ and then 
from u/ to Mj.* 

2 Prove that = 3, = » 

3 Evaluate the components of the sixth-order tensor 

^abe ®iin) - , ( f ' ^ ' 

4 Prove that ^ ^ 

5 Prove tliat " . 

^tka^mvs ~ ^stk^amp ^kst^pam 


6 Prove that if w,, v*, are vectors^ 

[W [u. v]]^ = u„ (VfW.) - v„ {u,w,). ^ 


Mm [». M’lm 

— ®«fcmM, 

A (m) = 

1 Wll 

Mu 

Mu 


1 'Uji 

Ma 

M» 


! 

Msa 

MjS 


7 If 


prove that 

*t}k ^ (**) ^ 

= 'hnn A (w), y 

6A («) = 

8 Use Ex 7 to prove the rule for the multiphcation of deter- 
mmants 

‘Imfl'lmn A («) A (v) = 6A (u) A (v) = 6A (uv), 
where (««),» = W|j‘’to- 



CHAPTER II 


\ 


GEOMETRICAL APPLICATIONS 

The displacement from any point to any other obviously 
constitutes a vector. The distance between the points is a 
scalar. If x ^ , are the coordinates of the points and r the 
distance between them, 

= (j/, - XiY (1), 

the square on the right indicating the scalar product of the 
vector into itself. Also the quantities (y, — a:<)/r constitute 
a vector. 

If we take a fixed point a, and consider points given by 

«. = «.+ hr (2), 

where r is a variable scalar and the are constants such 
that 

= 1 (3), 

{X, - a,Y = ( 4 ), 

so that r is the distance of x, from a, . If we take another 
point y, such that 

y, = a, + l,s (6), 

{y, - a,)^ = ( 6 ), 

{y. - = (« - ry (7), 

and therefore the distances between a,,x,, and j/^ are such 
that the sum of two of them is equal to the third. Thus the 
pomts are on a straight Ime, and (2) gives the equations 
of the hne in terms of the parameter r The h are the 
direction cosines of the hne. 

If we take two lines through given by 

Xt= ai + hr 

yf = af + mtS (f), 
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the distance between Xf and is given by 


17 

/ 


(!/i - = (-sw, - 

= s* + r® — 2rsZ,nii (10) 

But this quantity is also equal to s® + r® — 2rs cos 6, where 
6 18 the angle between the hnes Hence the angle between 
two intersecting hnes is given by 

cos 9 = IfW,, (11) 

If two hnes have the same direction cosines they are said 
to be parall^ If two hnes do not intersect we can take a 
line through any point on one of them parallel to the other , 
then this hne is inchned to the first at an angle given by ( 1 1 ) . 
We can then use (11) to determme a umque quantity 
associated with any two hnes, which we may caU their 
inclination, whether they intersect or not 

If we have a hne given by (2) and y, is a point outside it, 
the hne joining and y, subtends a nght angle at r, if 

iy. - «»)* = {x, - a.)® + {y, - 


= + (y, - a, - rl,y 

= (y. - a,)* - 2rl, (y. - a,) + 2rH.^ (1^), 

and therefore 

r = h {y, - a,) (13) 

This gives the projection of the displacement y, — a, on 
the hne. The foot of the perpendicular is 

cti + — a, + I, lie (Vk ~~ ^k) 

Evidently r in (13) will be the same for all pomts y, such 
that I, y, 18 constant Hence 

l,y^ = s (14) 

represents a plane perpendicular to the hne ^ 

If we take two intersecting hnes given by (8) and (9), we 


J CT 
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can find the equation of the plane containing them as 
follows. If this plane is 

»,z, = p (16), 

this equation must be satisfied by a;, and for all values of 
r and a. Hence 

n.at^p (16), 

nM = 0 ' (17), 

tiitni = 0 (18), 

and from (16) and (16), 

n, (z, - a.) = 0 (19) 

Then (17), (18), (19) are three homogeneous equations m 
the Jii , and can be consistent only if 

^tJcm (Zt ®t) ^k‘^m ~ ^ (20). 

This 18 the equation of the required plane Also the Uf are 
proportional to 

= ih‘^3 ~ ~ h'^> ~ ( 21 ) 

But n.* = 1 (22), 

(lj»i3 - + (Igmi - liin^Y + - h'frhY 

= (fi* + + is*) {m^ + + rug*) - (ij/ni + + is^s)* 

== 1 — cos* 0 

<= sin* 0 (23) 

Thus sin 0 n, = ± (24) 

The ambiguity in sign corresponds to a general one in 
specifying the parameter r of a point on a hne If r in (8) ^ 
IS taken negative, we get a pomt on the Ime on the opposite 
side of a, from those given by positive values of r But if 
we reverse both r and the Z, we still keep Z<* = 1, and we 
stdl have the same point We may take either direction 
along a hne to be that of r mcreasmg, if we reverse the 
direction the signs of all the Z, are reversed for the same 
pomt 
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For any point on (8), 

~ ^tkm (^k "i" ^k^) ^ ^ikm^k^m (^^)* 

This IS a constant vector for all values of r, and may there- 
fore be considered as a property of the Ime We denote it 
by 1/ Then we have six properties of the Ime given by 
li, 1/. These are coordinates of the line They are connected 
by two relations, 

= 1 (26), 

lilt = — 0 (27). 

The I,' have a geometrical mterpretation Thus if we con- 
sider the plane 

^2^3 ^3^2 ” (28), 

this plane passes through the hne. Also if x.^ = Oj,, — a^, 

(28) IS satisfied for all values of Xi, and therefore if we take 
a hne through a, parallel to the Xi axis, (28) represents the 
plane through (8) and this hne Two such planes determine 
the hne, and therefore the I, and f,' together determine the 
hne 

If we have two non-intersectmg hnes given by 

a;, = a, + rl, (29), 

y. = P. + sm, (30), 

the line y. = )3, + rl, (31) 

passes through and is parallel to (29) The plane mcludmg 
(30) and (31) is, by (20), 

( 2 . — ^t) ( 32 ) 

This therefore represents a plane through (30) parallel to 

(29) The plane through (29) parallel to (30) is 

e.km {z, - a,) hm„ = 0 (33) 

The distance between these planes is the projection of the 
hne joining any two points on them upon a hne perpen- 


2-2 
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mcular to both. H a line perpendicular to both has direc- 
tion cosines the shortest distance d between the hnes 
18 therefore given by 


d sin 0 = (^, — a,) n, sin 9 

= ± (34) 


But — ~ ‘^k t ^Ikm<^ik ~ (36)) 

and therefore 

d sin 0 = ± (— 

so that, apart from the ambiguity in sign, 

d sin 9 = l,m,' + w.Z/ (36). 

Thus the shortest distance is directly expressible in terms 
of the coordinates of the two hnes. 

Now consider two mtersectmg hnes 

x, = a, + rl,, y,-= a, + sm, (37) 

The area of the triangle formed by a, , x, , y, is ^rs sin 0 
The projections of these points on the plane x^~0 are 
(0, ffu, flg), (0, Xj, Xg), (0, ^ 2 ) Vs) u-nd form a triangle whose 
area is 



1 

Ua 

«3 

II 

«2 

Og 


1 

Xg 

Xg 

0 

l^} 

IgT 


1 

Vi 

2/3 

0 


nigS 


= Irs {Ignig - Igirig) 

= ^rsn-i sin 0 (38) 

Thus the projections of a triangle, and therefore of any 
plane area, on the coordmate planes are in the ratios of the 
direction cosines of the normal to the planes A plane area 
can therefore be treated as a vector whose components are 
proportional to the direction cosmes of the normal. 
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If we have a line given by 

x, = a, + rl, (39), 

and jSf IS a point not on the line, let us suppose the point j8, 
turned through an angle 6 about the line. The foot of the 
normal from j8, to the hne is given by 

= I, - a.) (40), 

and therefore the displacement from the foot of the normal 
to IS equal to 

A - (a. + l,h (A - «*)} = (A - a,) - hk (^t - «it) 

= k {k ~ U,) — k (^k ~ 

(41). 

The magmtude of this displacement, p, is given by 
p^ =L (^, _ a,)* _ {k ()8. - a.))2 

= k^ (^. - a,)* - kk (i8. - a,) - a,) 

= h {k W, - «,) - k (42), 

the i being needed because m the double summation each 
pair of values of the suffixes would occur twice 
The plane through and the hne is 
n, (z, - a,) = 0, 
subject to n,k ~ 9, 

n, iP, - a,) = 0 , 

and IS therefore 

e,km (z. - «.) k iPm - “m) = 0 (43), 

while the n, are proportional to €,^„,k iPm — «m) and there- 
fore equal to ± e.^k (Pm ~ a^)lP 

If now we turn p, through an angle 9, it receives a dis- 
placement p (1 — cos 6) along the normal to (39), and a 
displacement p sin 6 along the perpendicular to the plane 
(43) If it goes to y,, we have therefore 
y. - /S. = - (1 - cos 9) k {k iP, - «.) - h iPk - a*)} 

± sin ^ e,^k iPm - “m) (44). 
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If 0 IS a small angle and we neglect 6^, the displacement is 
simply 

Y, - Pi = ±6 - am) (46), 

and in particular 

Vl — = ± {^2^ (^3 ~ ~ (^2 ~ *^2)} (^®)> 


and from the additive form of this equation we see that the 
displacement is the sum of those given by separate small 
rotations 1,6 about axes through parallel to the co- 
ordmate axes Conversely, displacements due to small 
rotations about axes through a pomt can be added vecton- 
ally as if all were apphed to the system in its original 
position, and give the same total displacement as if they 
were compounded into a single rotation about an axis by 
the vector rule We still have, however, to estabhsh a sign 
convention. We decide that 6 is to be taken positive if a 
turn about the axis of is from towards Xg . This would 
make 


Yl~ ^ (p2 ®2) > Yz p2~ ^ {^1 ~ ®l) 

With — I 2 = 0, Ig = 1. Hence 

Yt ~ Wm ~ ®m) 


If we write 


l]c6 — UJic, 

Yi ~ ^tkm^k i^m ^m) 


(47) 

(48) . 


(49). 


For mstance, if = 0, we have 

Yl ~ ~ ^2^3 ~ ^3^2) Y 2 ~ Pz ~ ^3^1 ~ ^1^3> 

YZ ~ ^3 ~ ^1^2 ~ ^2 Pi (® 9 ). 

For fimte rotations we return to (44) or (41), keepmg the 
positive sign in the second term. We may transfer the 
ongin to a, to save wnting The coefiBcient of in y, — 

IS 


bik = (1 - cos 6) l,lk - sin 


( 61 ). 
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iox k ^ = i, the coefScient is — (1 — cos d) (1 — 

Thus 

6,* = (1 - cos 6) {l,h - - smd c.* (52), 

where is the antis 3 raimetncal tensor corresponding to 
l„, namely 

/ 0 h ~ h\ 

(-Is 0 lA (63). 

\ Is — li 0/ 

Thus the displacement is represented in general by , 
where 6^;^ is a tensor of the second order, expressed as the 
sum of symmetrical and antisymmetncal tensors The anti- 
sjmimetncal part is seen to be of the first order of magmtude 
in 9 and the symmetrical part of the second order. 


EXAMPLES 

1 Given that the general quadnc surface is 

S = + J5,a:, + C = 0, 

prove that the locus of the mid-points of parallel chords is a plane, 
and find the condition for this plane to be perpendicular to the 
chords 

2 Find the condition that the hne 

+ P = 0 

may touch the quadric 5. 
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PARTICLE DYNAMICS 

The essence of particle djmamics is that the second de- 
rivatives of the coordinates of a particle with regard to the 
time are equal to functions of the position and velocity of 
the particle with reference to neighbourmg particles The 
relations therefore provide a set of differential equations 
to determme the coordinates The equations for any 
particle can be put in the form 

wix, = SZ, (1), 

where m is the mass of the particle, Z< are the forces due 
to the other particles, and the summation is for all the 
other particles* It is a matter of experiment that this 
form holds when the axes are a certain type of Cartesian 
axes, which we call non-rotating and unaccelerated, or m 
brief dynamical 

If we take a different set of non-rotating axes at the 
same origin we have 

x,' = a„x, ( 2 ), 

and since the axes are not rotating the are constants 
Hence by differentiation 

Xj' = a,iX, (3), 

x/ = a„x. (4), 

and therefore velocity and acceleration are vectors 

The force Z, on the particle due to some other particle 
IS measured by the contribution to mx, due to the other 
particle, that is, the part of wx, that would disappear if 
the other particle was removed to an infinite distance It 
• Cf Jeffreys, Scienttfic Inference, Chapter vin, for a fuller analysis 
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is actually measured by the acceleration If we have the' 
three acceleration components of the particle 1 due to the 
particle 2, we can find the acceleration in the direction of 
Xj by the formula (4), and the form (I) will become 

mx; = SX/ (6), 

provided that we define X,' by the rule 

AY = a, A. (6) 

The meaning of the force in any direction requires de- 
fimtion in any case if we define it by its relation to the 
acceleration component in that direction, it follows auto- 
matically by (6) that force is a vector It follows that the 
sum of any number of forces, obtained by adding the 
components separately, is also a vector But the practical 
importance of the idea of force arises equally from the fact 
that in many cases the force components are known once 
for all from experience as functions of the coordinates and 
velocities The form (6) then says, as a general pnnciple, 
that the forces in any direction are additive 
If we form the contracted or scalar product of rax, and 
SA, by the vector a, , we get 

= 2A, r, (7), 

the left side of which is ^ (The square of course 

implies the product of x, by x, and therefore the summation 
for % — 1, 2, 3 ) We write 

T = = \m + x^^ + (8), 

and call T the kinetic energy of the particle Then by in- 
tegration with regard to the time from to we get 

\T =2 X,x,dt (9) 

But in any short interval of time di, x, dt is the increase 
of X,, namely dx. Hence the right side is equal to 
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S j Xtdx, taken from the initial to the final position of the 

particle We call the scalar X^dx, the work done on the 
particle hy the force X, m a small displacement dx^ . Then 

S I Jl , dx, IS the total work done by all the forces on the 

particle during the motion We have therefore the scalar 
relation that in any motion of a particle 

Increase of kinetic energy = Work done on the particle 

( 10 ) 

In a system of particles we may add upljthe equations (10) 
for all the particles We now take T, the kinetic energy of 
the whole system, as the sum of those of aU the particles. 
Then we get 

T = i:\mxA =22 X,dx, (11), 

L J<o L Jio Jh 

the double summation implying summation for all pairs of 
particles , the first summation is for the particles producing 
the forces Jl, and the second for the particles acted on It 
may happen that the X, are all functions of the coordinates 
alone, and not of their velocities, and that provided the 
initial and final values of the are the same the integral is 
the same however the x^ may vary in the interval If so, 
the mtegral is the difference between the values of a certain 
function U for the imtial and final positions of the system , 
we call the system conservative and V the work-fiinction 
Then (11) becomes 

or T — U = constant (13) 

This IS the equation of conservation of energy The quanti- 
ties on both sides are scalars — C7 is often denoted by V 
and called the potential energy 
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If U exists, then in all possible small displacements of 
the system 

SSZ.da;. = dU (14), 


where the dxt, 3n in number for n particles, are all in- 
dependent The coefficient of dx, for any one particle is 
SX,, the total force in the direction of a;, on that particle. 
Hence the force on a particle in the direction of is 
dUjdx,, which is of course a vector 

We can also take the scalar product of (1) by any set of 
small quantities Sa:, whatever that constitute a vector. 
Then 


ma:,8a;, = SX,Sx^ (16) 


IS a scalar relation, but as the 8a;, are arbitrary we can 
equate their coefficients and regenerate the equation (1) 
If we now mtegrate with regard to t from tg to ij we get 

f mx,Sx,df=[ EX,SXfd( (16) 

Jh Ju 

The left side, on integration by parts, gives 

mx^hxA — mx^ j hx,dt (17) 

Jfo Jh 

But we can consider the x, -f 8x, as coordinates of a par- 
ticle m a motion diflenng shghtly from the actual one , that 
18 , at a given value of ^, the coordinates are x, -t- 8x, instead 
of X, Then 

a -a*' 


= X, -h 8x, — X, 

= 8x. (18), 

since X, -I- 8x< is simply the varied velocity or rate of change 
of the varied coordinate x, -(- 8x,. Then 


f mx. Sx, dt = 

Jto dt 


f 

JU 

!'■ 


mx,8x,dt 

8 (\mx,^) dt + 0 (Sx,)^ (19) 
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Then for every particle, to the first order, 


{S ( Jmac,*) + SX,Sa;,} dt = 

to 



( 20 ) 


We can now add up these equations for all the particles of 
the system If U exists we can express the result in the 
form 


S (T + t/) = S 

Ui 


OTX.Sa;, 


( 21 ), 


or, since the hmits are not varied, we can move the 8 out- 
side the integral If the vanation is such that the imtial 
and final positions of the system are unaltered, 8z, = 0 
when t = or ti, and we have 


8 \^\t+ U)dt=0 
iu 

to the first order in the vanations of the path This is 
Hamilton’s principle 


EXAMPLES 


1 If the Cartesian coordinates of every particle of a system are 
known functions of a set of generalized coordmates g, , prove that 




Hr’ 


where the summation convention is also imderstood on the right 
side Deduce Lagrange’s equations for the case where the are 
all independent 


2 The equations of motion of a particle are widJ, = X, — 
where k is constant Prove that 

d 

2T= -x,X, + 

Hence show that for a system m periodic motion, or in one slowly 
changing its state, on an average over a long time, 2r = — 2a;, X, 
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DYNAMICS OF RIGID BODIES 

A ngid body is one such that whenever it is displaced the 
distance between any two particles of it is unaltered. Since 
three particles A, B,C are in a straight hne if the sum of two 
of the distances AB, BC, CA is equal to the third, it follows 
that straight lines are unaltered by displacements of a ngid 
body Smce when A,B,C are not in a straight hne the angles 
of the triangle ABC are deterimnate functions of the three 
sides, it follows that all angles are unaltered by displace- 
ments of a rigid body If three hnes meeting at a point and 
fixed in the body are mutually perpendicular before dis- 
placement, they are still perpendicular after the displace- 
ment 

The equations of dynamics in the form mx, == X, are 
true with respect to dynamical (that is, non-rotating and 
unaccelerated) axes Let any particle Q of a body have 
coordinates x, with reference to dynamical axes at 0 Then 
let the body be displaced in any way, and let the particle 
have the new coordmates x/ We require the relation 
between x,' and x, Consider a particle of the body, P say, 
whose coordinates before and after the displacement are 
a, and a,' Put 

x, = a, + y„ x/ = a,' -t- y,' (1) 

Then y, and y/ are the coordinates of Q with respect to 
axes at P parallel to the dynamical axes before and after 
the displacement Also if we imagine the original axes at 
P to be specified by the particles on them, these particles 
in the new position still specify a set of rectangular axes, 
with respect to which the coordinates of Q are still y, If 
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the cosine of the angle between the axis in its new posi- 
tion and the yt axis in its old one is o,* , we have therefore 

y,' = aocVk (2), 

whence x^' = a/ 4- {x,, — Oj.) (3) 

The displacement of Q is 

o/ - at + a,^y^ - y, (4) 

We can prove that there is a straight hne of particles m the 
body such that the displacement of Q is the same as that 
of P. For such a particle we should have 

aaVk - y. = {a,k - Kk) 2/i = 0 (5), 

and these three homogeneous equations in y* are consistent 
provided 

= 0 ( 6 ), 

that IS, a^i — 1 0,2 0,3 =0 (7) 

0 „l ^22 ^ Og.q 

®31 ®32 ®33 ^ 

Now the determinant 1 1 0 ,,^ | j is umty, and each element in 
it IS equal to its first minor If we expand (7) we get 

II II “ {(®11®22 ~ ®12®2l) + (®22®33 ~ ®23®^32) 

+ («33«u - aaiOis)} + («11 + ^22 4- 033) - 1 = 0 (8), 

Since the terms all cancel Hence (5) have an infinite 
number of solutions, all proportional, and the points there- 
fore he on a straight line 

If we take any two planes through this hne, the angle 
between them is the same after displacement as before, 
and therefore all planes are rotated through the same angle 
Thus any displacement of a rigid body is equivalent to a 
displacement of a particle of it combined with a rotation 
about a hne through that particle If the angle of rotation 
IS 6, we have, by companng (5) with (52) of Chapter ri, 

^tk = ~ ^ik 

= {hk - S,t) (1 - cos 0) - sm 0 (9), 
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where Z* are the direction cosines of the axis of rotation 
Thus the are determined in terms of this axis and the 
angle of rotation. Conversely, if we take the symmetrical 
and antisymmetncal parts, 

~ ^tk) (1 ~ COS 6) = ^ (b/i + b/ci) (10), 

i (etim - e,km) Zm sm 0 J (6,* - 6^,) (11) 

But for a given m, with ikm m cychc order, 

^kim ^ikm ~ ^ 

Thus Z„ sm = - I (6.fc - 6^) 

and 

sin* (9 = J {{bi2 - b^i? + {b^s - 63^)* + {b^i - 633)*} (14) 

Now suppose that the displacement is small Then 
have nearly their values for zero displacement, that is. 
The direction cosines of the y* and axes in their new 
positions are a,* and a,i, and thus, d k ¥= I, 

a.fca., = 0 (15), 

and, if k = I, =1 (16) 

But in (15) for t — k, a,* is nearly 1 and a,i small, for t — I, 
a,i 18 nearly 1 and a,* small, for t not equal to k or Z, both 
o,j. and o,[ are small Hence to the first order 

aj* + au=0 (17). 

If ^ = Z, then for i¥^kov I both terms of (16) are small of the 
second order, and therefore for i = k or Z, a,i = 1 + a 
second order quantity Thus to the first order 

^tk “ ^tk (16), 

where 6*4 is an antis5mimetncal tensor 
The displacement of Q is 8a, + a.^yj, — y,, where Sa, is 
the displacement of P, and is a first order small quantity. 


(12) 

(13), 
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Now let X, be the total force acting on the particle at Q. 
In the displacement the work it does is the scalar 

SIF = X, (8ai + a^yf. - y,) 

— X, (Sa, — bficyic) 

~ X, Sffl, (X,yic + Xj.1/,) (X, 2/j. Xjji/,) 

(19) 

When t and k are interchanged, and X.i/t — X^y, are 
reversed in sign, while X^y^. + X^y, is unaltered. Hence 

SW = X.Sa, + J (i/,Xi - y„X,) b,„ (20) 

Now wnte y^X^ - y^X, = L,„ (21), 

so that Ljic IS an antisymmetncal tensor The second term 
in (20) IS the sum of mne terms, of which three are always 
zero and the others equal in pairs If we replace 6,*. and 
by the associated vectors, we have 

SW = X,8a, + L^b^ (22), 

where the Sa, and b^ are all independent of one another 
and the same for all particles of the body Hence if we add 
for all particles of the body 

SW= (i:X.)Sa.+ (2L„,)6„, (23), 

and the total work done in any given small displacement is 
determinate if we sum up the forces acting by the six 
expressions 2X,, Further, the contributions to 

SX, , from the internal reactions are zero This follows 

at once if these reactions consist of equal and opposite 
forces between pairs of particles along the Ime joimng 
them, and also has the justification that it leads to correct 
results Then we may restrict SX,, to the contribu- 
tions from the external forces This is d'Alembert’s Prin- 
ciple 

In the hmiting case of continuous motion, we may con- 
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sider the displacements that tak^ place m a short time 
interval St, and put in the hmit 

Sai = u,St, (24), 

and call the velocity of P and a>,* the angular velocity 
of the body. Then the velocity of Q is 


*. = M. - a).iy4 (25). 

Now consider the centroid O, with coordinates de- 
fined by 

(2m) = 2ma:, (26). 


It IS usually assumed without proof that O is fixed in the 
body, though this is not obvious But suppose that the 
particle at G with coordmates a, m the ongmal position of 
the body goes to O' when the body is displaced, its new 
coordinates are a/, given by 


j-i — a, + o,* — cik) (27), 


and the coordinates of the new centroid G" are a , ", given by 
(2m) a," = 2ma;,' (28). 


We have to show that G' and O" coincide. We have 


(2m) a," = 2m {a,' + (x^ - o*)} (29), 

and therefore 


(2m) (a,* - a,') = 2m {a,' + (x* - %)} 

- (2m) (a.' + o.fc (Jfc - o*)} 

2ma,]^Xj^ “ 2mu,j^aj^ 

= 0 (30). 

Thus the particle ongmally at Q is displaced to the new 
position of the centroid, and therefore the centroid is fixed 
m the body 

Now return to the equations of motion of a particle of 
the body 

mx, = X, (31) 


J C T 


3 
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By addition we form from these the equations 

"Lmi:, = SX, (32), 

and by cross-multiphcation followed by addition for all 
the particles of the body, 

Sm {XiXk - x^x,) = 2 (x.Zfc - a;*X,) (33). 

On the right sides of (32) and (33) the contributions from 
the internal reactions are zero, by d’Alembert’s prmciple. 
Also 

Sma;, = Sma;, = (Sm) a, = ilf?, (34), 

where M = Sm (35), 

the total mass of the body Also if 

/i.j. = 2m {x,Xk - XkX,) (36), 

we can reduce (32) and (33) to 

MXt = 2X< (37), 

~ A., = 2A., (38) 

These are the fundamental equations of rigid dynamics. 
The three mdependent A,* are expressible in terms of the 
three mdependent co.jt by (36), and we have therefore six 
differential equations for the a, and cu.^ , which determine 
them m terms of the imtial conditions if the external forces 
are given. The motion of the body is therefore determinate 
The principle of virtual work follows immediately For 
if a body is initially at rest, a, and A,* are zero, and the 
condition that they may remain zero is that 2X, and 2^,^ 
shall vamsh. But this imphes, by (23), that in any small 
displacement of the body the work done by the external 
forces IS a small quantity of the second order in the dis- 
placements. Conversely, if there are six independent 
possible small displacements such that the external forces 
do no work, to the first order, m any of them, the co- 
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efficients of the independent 8a, and of (23) must all 
vamsh. But these coefficients are the and Sim* 
therefore the a;, and A,* do not vary with the tune. The 
vanishing of the work done by the external forces in six 
independent small displacements is therefore a necessary 
and sufficient condition for equihbnum 

The equations (38) may be put m another form. If we 
consider any moving ongm O', not necessarily fixed in the 
body, with coordinates o, , we can write for the coordinates 
of a particle with respect to O', 

y, = a:, - a, (39) 

If any vector associated with the particle, such as its 
velocity, momentum, or acceleration, or the force acting 
on it, has components u, , we may form the antisymmetncal 
tensor y,u^ — and call it the moment of the vector 
about O' From (32) and (33) we can form the equations 

Sw (x,Xk — x^x,) — {a,'LmXi. — a*Swa:,) 

= 2 {x,X, - x,X,) - (0.2X, - a,2Z.) (40), 

where only the external forces make any contribution to 
the right side. But by (39) this is equivalent to 

(j/<** - ykX.) = 2 (y,X^ - yicX,) (41). 

Therefore the sum of the moments of the mass-accelera- 
tion products about any ongm is equal to the sum of the 
moments of the external forces about that origin. 

If we denote the moment of momentum, or angular 
momentum, about O' by A,*.', we have 

Kk = Sw (y,Xi - y^x,) (42), 

and 

^ A,t' = 2m (y.Xj - y^x,) + 2m (y,x^ - y^x,) 

( 43 ). 


3 -* 
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The second term may be wntten 

Sm {(Xi — cf<) x^ — (Xi — a*) i,} = — 2m (a^x^ — a^x^) 

= - if (a/Xf: - a^a,) (44) 

Thus (41) IS equivalent to 

^ + if (Uja,. — tt^a,) = S (i/{X)c — Vic^i) (46). 

In many important cases the second tensor on the left 
vamshes identically This is clearly true if O' is fixed, when 
the time derivatives of the o, are zero, when the centroid 
IS fixed, and, if both are moving, if the vectors d, and a, 
are parallel, that is, if the velocity of the moving origin O' 
18 parallel to that of the centroid. The most important case 
IS where the moving origin is identical with the centroid, 
when the last condition holds automatically These terms 
also disappear if the moving ongin is an instantaneous 
centre of rotation always at the same distance from the 
centroid They vamsh for a sphere or circular cyhnder 
rolhng down an inclined plane, but not for a rolling eUiptic 
cyhnder 

If the movmg origin is the centroid, 

Kk = Sm {y, (.fjt + y^) - y^ (x, + y,)} 

= Sw (y.yft - Vky,} (46), 

since 2my, = 2my^ = 0, 

by the defimtion of the centroid Thus the angular mo- 
menta about the centroid are expressed completely in 
terms of positions and motions relative to the centroid, and 
the formulae for them have the same form as those for the 
total angular momenta with reference to a fixed origin 

By (25) y, = - Ok.mVm (47), 

m being here a dummy suffix Substituting in (46) we have 
Kk = - Sm iu}ic„y,y„, - (Oir^ykyJ) 


( 48 ), 
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where is the symmetrical tensor defined by 

(49), 

and depends only on the masses and positions of the par- 
ticles of the body In the three-dimensional case, with i 
and k unequal, m must be equal m turn to z, k, and the 
other value. Then, for instance, 

Aj 2 = (-^ai^ll “t" -^22^12 “b -^23^13) 

— -h 5 j 2 a)j 2 + B13OJ23) 

~ (-®22 "I" -^ll) “*12 •®13“^23 — -®23“'31 (®®)- 

If we replace the antisymmetncal tensors by the associated 
vectors this takes the form 

A3 = (Bii -f- B32) 0J3 B13W1 B23W3 

. — (.Bu + ^22 -I- .B33) OJ3 — (B13CU1 + B33OJ3 -I- .6330)3) ( 51 ), 
and in general 

K' = ( 52 ), 

where . 4 ,^ is the sypimetncal tensor given by 

= (Smyi*) S,„ - I,my,y„, ( 53 ). 

It evidently corresponds to the system of moments and 
products of inertia given in ordinary treatises on dynamics , 
^11, .422, -^33 *^re the ordinary moments of inertia A, B, C, 
but A23, .dji, .4 i 2 are equal and opposite to the ordmary 
products of inertia F, G, H. 

It should be noticed that this reduction is characteristic 
of three dimensions; in a higher number of dimensions 
there is no analogous simphfication of the form ( 48 ) 

The equations of motion of a rigid body then take the 
form 


Mt, = SZ. 

= or ^A„' = SL„,' 


( 54 ) , 

( 55 ) , 
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where is the moment of the external forces about the 
centroid 

It may happen that one point of the body is fixed. In 
that case we may take the origin at that point. Then in 
addition to the known external forces there is a reaction 
at the origin which can if reqmred be found from (37) , but 
the reaction has no moment about the origin and the 


motion IS given by (38). But m this case 

~ ^tm^m (®®)» 

and we find, by a process analogous to the last one, 

where Z),„ = YimXiXm (58) 

In terms of the associated vector, 

(59), 

where (7,„ = (Smx^^) 8,^ — 'Lmx,Xm (60). 


The G,n correspond to the moments and products of inertia 
about the origin They can be expressed in terms of those 
about the centroid, for 

= {Sm (Jt + VkY) - 2m (J, + y,) + y„) 

= (2m.Ts2) 8,„ - 2m.T.x„ 

+ ( 61 ), 

the terms Imear in the y’s vanishing by the definition of the 
centroid 

The relevant equations of motion then take the form 

^^h„, = 2L,„ (62), 

where the are the moments of the external forces about 
the origin, the reaction at the ongm making no contri- 
bution. 
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The kinetic energy of the body is given by 
2T = Smar,* 

= Sm (i. - 

= Swi.t,2 + 

= ito,„ (63) 

In three dimensions co.^ajj^ is zero unless i is different from 
both k and m If i = 1, i = 2 , m = 3, or if i = 1, jfc = 3, 

Wl ^ 2) Ct)22^13 — ^2^3* ^ ”” 1 j A/ 2j ^12^ ^ ^3^> 

and has coefficient Thus enters with a coefficient 
■®22 + -® 11 J or ^ 33 . But <i> 2<^3 coefficient — 2 B 23 or 

^23 + -432. Thus in all 

2T = (64) 

When a point of the body is kept fixed, 

2T = C)cm*^K^m ( 65 ) 

We notice that the linear and angular momenta take the 
forms 

JUt, = ^ constant), 

CcCj (TCOj 

ST 

for origin at a point of the body held fixed 

In one respect the forms (55) and (62) are inconvement 
They involve the tensors A,„ and which depend on 
the and y, and therefore in general change as the body 
rotates It is more convenient to use such axes that the 
relevant tensor in the particular problem is constant To 
achieve this the axes must rotate, and then are no longer 
dynamical axes Suppose then that we have a set of axes 
x/ rotating in any way, and that their direct iqp cosines 
with respect to the dynamical axes a:, are a,j All the usual 
tensor relations hold for transformations from the x,' 
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system to the x, system and mce versa; but the while 
specified for each instant, are now functions of the time 
If M, IS any vector, we have 

Mj' = Ut = at)Uj (67) 


Then 


dui 

dt 


IS another vector; and 

dut _ dUi , da„ 
dt dt ^ ^ dt 


(68) 


The component of this vector m the a;,' direction is 


du/ 


+ 0,1% 



(69) 


The first term of this is 


du/ _ 
^^‘-dt- 


du/ 

dt 


(70) 


da. 


Also 18 the X, velocity of a point at umt distance along 


the x/ axis, that is, a point with coordinates o,, (t = 1, 2, 3) 
with reference to the dynamical axes This velocity is 
— dj^a^j, where is the antisymmetncal tensor expressing 
the rotation of the ngid frame consisting of the moving 
axes. Then 


0,!% 


, dUf, 

dt 




= - B,/Ui 


and the required component is 


dt 


- St/u/ 


(71) , 

(72) 


If we use instead of 9j/ the associated vector, the three 
components become 

(%' - u/d/ + u/ 9/ , u/ — u/ 9/ + u/ 9/, 

- u/e/ + u/9/) (73). 
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If for instance the components of a displacement along the 
Xj axes are . the formula (73) gives the components of 
the velocity along these axes , if Uj are the components of 
velocity, (73) gives the components of acceleration, if «/ 
are the components of angular momentum, (73) gives the 
components of their rates of change with reference to 
dynanucal axes, and these are equal to the components of 
the moments of the forces actmg 

EXAMPLE 


Prove that 
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EQUIVALENCE OF SYSTEMS 
OF FORCES 

An external force Xt acting at a point Xi of a rigid body 
produces d 3 manncal effects summed up m the vector Xt 
and the antisymmetncal tensor a;, A* — x^Xi A force has 
a hne of action , that is, if we take its resultant R given by 

( 1 ) 

we can define a direction Z, by 

A. = Rh (2) 

By convention R is always taken positive. If a force A, 
acts at the point x, + , where r is arbitrary, we have 

(x, + rZ.) Xy. - (Xfc + rh) A, = x, A* - x^X, (3). 

Thus the dynamical effects are the same if the force A, 
acts at any point of a hne through x, with direction cosines 
proportional to A, This hne is called the line of action of 
the force, and the force can be said to act along it 

If a force has magmtude R and acts at x, m the direction 
Z, , we have 

A, = Rl, (4), 

Rik — R ^kl%) (5)> 

or = Rl„' (6), 

in the notation of the coordinates of a hne Thus the A, 
and are the products of the resultant into the six 
coordinates of the hne of action 

The moment of the force A, about a point a, is 

(x^ — a,) Xj^ — {xj^ — a*) A< = R — o,Zj. + o*Z,) (7). 
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If we take the associated vector, its component in the 
direction is 

■R - e,kmin^ath) = -R (8), 

where m,' are the other coordinates of the hne through a, 
in the direction to, . We notice that 

mj,' + l^m,' = dsmd (9), 

where d is the angle between the two hues and d the length 
of their common perpendicular We may call Rd sin 6 the 
moment of the force about the line (mt, to,'). 

By d’Alembert’s principle, the motion of a body is un- 
altered if to the forces acting on it we add two equal and 
opposite forces acting at the same point, or, by (3), along 
the same hne of action. 

Now consider a pair of equal and opposite forces X, and 
— Xt acting at points a, and b. They clearly make no 
contribution to SX, Their contnbution to is 

(a^Xjc — UftX,) — (b^Xjc — 

= (a, - b.) X^ - [a^ - X< (10) 

Since a, and 6, are equally affected by any motion of the 
origin, the contribution of such a pair of forces apphed to 
defimte particles to both 2X, and is independent of 
the position of the origin Such a pair is called a couple, 
and its contribution to is called the moment of the 
couple 

If the vectors a, — 6, and X, are both perpendicular to 
a line with direction cosines n, , the components of are 
the components of a vector along this hne , while the forces 
act in the same plane perpendicular to the hne This plane 
IS called the plane of the couple, and the hne an axis of the 
couple Evidently equal couples in parallel planes are 
equivalent The magnitude of this vector is Rd, where 
R^ = Xf and d is the perpendicular distance of o, from a 
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line through 6, parallel to -X,, and its components can 
accordingly be written Rdn„ 

Any system of forces is eqmvalent to a force at an 
arbitrary pomt a, together with a couple. For with each 
A, actmg at a;, we can associate a pair of forces ± at . 
Then the system is equivalent to at a, together with 
a set of couples whose total moment is 

= S {(x. - a.) X;i, - (x* - o*) X,} (11). 

But this IS equivalent to a single couple ; for we have only 
to make proportional to the M„ and Rd equal to their 
resultant. 

If L., are the moments about the origin, 

^ik = 2 (x.Xj. — XjcX,) (12), 

and if for brevity we replace SX, by simply X, , 

M\k — ~ ~ (13), 

or i/„ = 2/„ - ea„a,X^ (14) 

Evidently X,“ is a scalar and independent of a,. Also 

Xrn^m ~ ^tkm^iXj^Xyn (i^) 

The first term is a scalar and independent of . The second 
is identically zero Hence X,^ and X^M( are scalar in- 
variants. 

The system is equivalent to a force X, and a couple M, 
at a^ . These vectors are parallel if 

Lm — ^ikm^t^k = pXjn (16), 

where p is a scalar length called the pitch These give three 
linear relations between the three a, and the pitch, and we 
therefore expect a single infimty of solutions. But if we 
take the scalar product of (16) by X„ we have 

Xm^m = pX„* 


(17). 
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SO that p 18 the same for all admissible values of . With 
this value of p, (16) represents three planes with a line in 
common. If we change a, to a, + aXi, where tr is a scalar, 
the left side of (16) is increased by 

~ ^tkm<^X^Xjc — 0 ( 18 )> 

and therefore if a, is one point on the line, all points on the 
line through u, parallel to satisfy the conditions This 
line is the central axis of the system , the system is equiva- 
lent to a force R along the central axis and a couple O 
about it This expresses the system as a wrench If we take 
the central axis as one of the coordinate axes, we have, 
since Xt^ and X^M^ are scalars, 

R^ = (19), 

GR = (20), 

and therefore G, R, and p are determined 

The system can also be reduced to a couple parallel to 
a preassigned plane together with a force For if S is the 
couple, and are the direction cosines of the normal to the 
plane, and if the force acts through a, , we have, for the 
moments about a, , 

^■m ^ikm^xXk (^1)) 

and also = Sn„ (22) 

We have three equations to determine the a, and 8. Agam 
there are a single infimty of solutions But if we take the 
scalar product by X ^ , we get 

RmXff^ €fj^jf^a^Xj^Xj^ 

= GR (23), 

so that 8 IS determined provided n^X^ is not zero, that is, 
provided the resultant force is not parallel to the plane. 
Then the equations 

^tkm^iXjg = 8nJJ^ 


(24) 
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determine a line parallel to JC, , which is the line of action 
of the force 

If If are the direction cosines of a hne through a<, the 
moment of the system about this hne is 

~ ( 25 ) 

= ImLm + h' (26), 

where Z*' are the other coordinates of the hne. If this 
moment vanishes the Ime is called a null hne of the system. 
If 6, IS another point on it, 

Lm (bm — Om) = ^,km<^i^k (bm ~ ®m) (27), 

which shows that 6, hes in a definite plane through o, . All 
null hnes through a point therefore he in one plane This 
plane is called the null plane of the point 

All null hnes in a plane pass through a point For the 
system can in general be reduced to a couple in the plane 
and a force whose hne of action mtersects the plane in one 
point Then the system has no moment about any line in 
the plane through this point, which is the null po%nt of the 
plane 

Any system is equivalent to two forces, one of which can 
be made to act along a given hne For let the lines of 
action pass through a, and 6, , and have direction cosines Z, 
and m, , and let the magmtudes of the forces be 8 and T 
Then we have six equations, 

X. = 81, + Tm, (28), 

~ ^tkm “!■ b^TjUi^ 

= 8l„: + Tm„' (29). 

The coordinates of the first line being given, these are six 
hnear equations to determine the six coordinates of the 
second hne and 8 and T But we have also 

= 1 , m,m’ = 0 


(30). 
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so that we have m general just enough equations Two such 
lines are conjugate lines Clearly any hne intersecting two 
conjugate hues is a null hne 
It can be shown easily that 

/S {XJi + L,l,) = L^Xi (31), 

= X,^ ~ 2Sl,X, + (32), 

whence (28) and (29) determine the coordinates of the 
second hne exphcitly. 


EXAMPLES 

1 A system of forces is reduced to a force at P together with 
a couple, P IS chosen so that the couple is parallel to a given 
plane Show that the locus of P is a straight line parallel to the 
central axis 

2. Show that 



IS a point on the central axis 

3. Two systems of forces are given by {Xf , L,), ( F, , K,) Show 
that X, Kf + F, L, is invariant 
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CONTINUOUS SYSTEMS 


In problems mvolvmg volume and surface mtegrals we find 
it convement to denote elements of volume and surface by 
dr and dS respectively, both are always taken positive. 
They are of course scalars The direction cosines of the 
normal to an element of surface are usually denoted by I, , 
in most cases the normal is drawn outwards from the region 
under consideration 

Green’s Lemma takes the form 



1 l,UidS 

(1). 

and we have the corollary, if 



u, = dVfdXi 

where T is a scalar. 

(2). 



(3), 

, r7, 9* 9* 5“ 

Where = g + + 

(4), 


and djdn denotes differentiation along the outward normal 
Stokes’s Theorem takes the form 


L 


u,dx. 



^^ 3 / \v^3 


= Jj dujdx^ds 


dus\ 

dxj 

+ h{^ 


dUg 

dxi 



dS 


(5) 


The integral on the left is round a closed contour C. On 
the right I, is the direction cosine of the normal to any 
element dS of a surface S that fills up the contour; the 
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integral is over S , and the sense of the normal is such that 
if the contour is described m the positive sense about any 
axis a;, , Z* is taken positive when the normal is m the direc- 
tion of X, mcreasmg. 

The gravitational potential of a distribution of matter is 
given by 

F = /S^ (6), 

where m is an element of mass, r is its distance from the 
pomt where V is to be found, and / is a constant equal to 
6 66 X 10 ~® when m and r are measured in grams and centi- 
metres and the time m seconds When the mass is distributed 
contmuously over a surface or through a volume, m must be 
replaced by adS or pdr respectively, where a and p are called 
the surface and volume densities The work in displacing 
a mass m' through a smaU distance is m'dV. 

The electrostatic potential of a system of point charges 
IS given by 

V = f^l ( 7 ), 

where e is a typical element of charge If e is m electro- 
static units of charge and r m centimetres, /is + 1 The work 
in a small displacement of a charge e' m the field is — e'dV 

The usual relations follow, that in free space in both 
cases 

= 0 ( 8 ), 

and m space occupied by matter of fimte density 

V*F = — 477/0 ( 9 ). 

Also we have Gauss’s Theorem 

11 dS = — 477/E'm or — 477/S'e ( 10 ), 

where the summation is over all the masses or charges 
withm the closed surface S. In crossing a surface where 

JOT 


4 
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there is a finite surface density tr, dVjdn has a finite dis- 
continuity — infa. 

All these relations are scalar in form. For a gravitational 
field the force on a small particle of mass m' is 

m'Xt = m'dV jdxf (11), 

while for an electrostatic field the force on a small charge 
€>* IS 

e'X, = - e'dVjdx, (12) 

The vectors X, are called the intensities of gravitational 
and electric force respectively The analogy m form be- 
tween (6) and (7) is constructed for mathematical con- 
vemence; the difference in sign in (11) and (12) embodies 
the physical difference that whereas two positive masses 
attract, two positive charges repel 

In a gravitatmg system we may construct a work- 
function 

(13), 

where the summation is over all pairs of particles m^, m^, 
and rj„ is the distance between nij, and m ^ The case where 
p = q IB excluded Then the force on the pth particle is 
d W /(Sa;,) j, The potential at m^, due to the other particles is 
Fj, = the accent indicating that the case where 

q = p IB excluded from the summation , and the function 
Fj, = 2W, since each pair of particles is counted twice 
m this double summation The function Vj, therefore 

plays the part of a work-function Similarly in electro- 
statics the function JSe^Fp plays the part of the work- 
function with its sign reversed, that is, of a potential energy 
These results may be generahzed to the case of continuous 
distributions ; thus we can replace these functions by 

W = ijjj pVdT+ i jjaVdS 


(14) 
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where the first integral is through all regions where p is 
finite and the second over all surfaces where a is finite , the 
first integral may therefore be taken through all space, 
excluding the surfaces where there is a surface density. It 
can be shown that with proper precautions about the de- 
finition of V the restriction p ^ q gives no trouble provided 
V 18 everywhere fimte. 

Now if we consider the mtegral through all space except 
thin lanunae surroundmg the surfaces where a is fimte, 





where the dn in the first integral is out from the region of 
integration and therefore towards the surface where a is 
fimte On the two sides of such a surface the values of V 
differ by an mdefimtely small amount, and for the two 
dV 

sides together ^dS — 4^fadS, by Gauss’s Theorem 

Hence the first integral in (16) is equal to — ^ jjvVdS 

taken over the surface and cancels the second term in (14) 
Hence 



through all space In consequence of this form we may say 
that the gravitational work-function, or the electrostatic 
potential energy, is jR^/Stt/ per umt volume, where B is the 
resultant of the appropriate intensity vector 

When the properties of the medium vary from place to 
place, V IS no longer of the form ^fefr, and V^F is no longer 
zero But a potential still exists, if a small charge e' is 


4-2 
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moved from a point P to a pomt Q, the work done is still 
the same whatever the route, and may be denoted by 
e' (Vj>— Vq). The treatment is suggested by the fact that 
two similar condensers with the plates at the same po- 
tentials, but with air between the plates of one and another 
material between those of the other, have charges m a ratio 
K dependmg only on the media dVIdn on the outside of 
the condenser bemg small compared with its value be- 
tween the plates, we infer that the charge per umt area for 
the same distribution of V is related to the discontmmty in 
KdVjdn, where K depends on the material This suggests 
in turn that Gauss’s Theorem must be replaced by 

\\K^d8^-infi:e ( 17 ), 

where the summation is for all charges inside S Then 
applymg this to the two sides of any surface we have 

and applying it to a region with a fimte volume density we 
have 

( 19 ), 

and as this must hold for all such regions. 

These equations are all homogeneous in V, p, a, hence 
the potential due to any set of charges is proportional to 
the charges if all are altered m the same ratio. Usmg this 
prmciple we can show by the usual method that the energy 
of a distribution is 

^SFe= ^|j|pydT+ ^ jjaPdiS 


(21) 
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The first integral, apphed to all space except thin regions to 
cut off surfaces where there are surface densities, gives 

The surface integral cancels the integral ^jjaVdS, and 

so that the energy can be considered equal to KR^j^nf per 
unit volume 

Magnetism may be treated similarly, starting from the 
assumption of volume and surface distributions of mag- 
netic pole strength, subject to the condition that the total 
pole strength m any sohd is zero, or we may regard the 
ultimate magnetic unit as the doublet, which explains the 
need for the restriction mvolved in the former method of 
treatment The potential at x, due to a doublet of strength 
M at the origm with its axis m the direction Aj is 

F = yMX,x,jr^ (24), 

and if the doublet strength per umt volume in a sohd is / 
in the direction A, we can mtroduce the intensity of mag- 
netization at the vector /A,, and say that the 

potential at a; < is 

F = yjJ|^.^3-t*dT (25), 

where Ai corresponds to the pomt and dr — 
y is a constant. The magnetic force m free space is 

dV 
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We may write 
F = 



dr 


= V 

= y 



IdAj 
r Sit' 


dr 



dA,i, 
Si. r 


(27) 


The potential is therefore eqmvalent to that due to a dis- 
tribution of magnetic poles l,A, per umt area over the 
boundary and — dA.jdit per umt volume through the in- 
terior. 

Within a solid special treatment is needed To define F 
or the force at x , , when x. is within a sohd, we must imagme 
a small cavity made about x. , the intensity of magnetiza- 
tion everywhere remainmg as before, and consider F and 
a. within it ; then the values of F and a. at x, are defined to 
be the hmits of those in the cavity when the dimensions of 
the cavity become mdefimtely small This process leads to 
little difficulty in gravitational and electrostatic problems, 
but in magnetism the hmit of the force is found to depend 
on the shape and orientation of the cavity The force in 
the cavity can be written 


^. = - 


dV 

dx. 


(28), 


where F is given by (27) , in the first mtegral the normal is 
inwards towards the cavity The contributions to from 
the volume integral and the outer boundary are of the same 
form as for gravitation, and give no trouble If the cavity 
IS a cylmder m the direction of the mtensity of magnetiza- 
tion, l.A. IS zero over the sides and equal to I, the resultant 
intensity of magnetization, on the ends Such a surface 
density over the ends m the hmit contributes nothing to F ; 
if the radius of the cylmder is small compared to its length 
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it also contributes nothing to X^, but if the cyhnder is 
of disc-hke form it contributes to X, Hence if we 

take F for the complete body and define a, by 


a. = 


^F 

3a;, 


(29), 


a, is the value taken by X, m a thm cylindrical cavity 
parallel to the intensity of magnetization The force in a 
flat cylindrical cavity with its generators in this direction is 

a, = «, + 47ryA,7 = a, + 47Ty-i4, (30). 

Evidently a, and a, are both vectors, the former is called 
the magnetic force and the latter the magnetic mduction 
The theory of susceptibility and permeability may then be 
developed as usual Also (27) shows that F is continuous 
across a boundary, but I, a, = dV/dn has a discontmuity 
— whence l,a^ is continuous across a boundary 

The mutual potential energy of two doublets M and M' 
at X, and a;/, onented in directions A, and A/, is 


IF 


K'M' 


dx,' 


( 31 ), 


where F is the potential at x,' due to the magnet at x , , this 
gives 


IF = yMM’K' (a. 



f ^ k ~ 

|73 


= yMM' 3 


yMM' n o D'\ 

■ ' g (cos € — 3 cos 9 cos 9 ) 


(32) 

(33) , 


where € is the angle between the axes of the magnets and 
6 and 9' are the angles made by the axes with the hne 
joinmg the centres. 
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If the second magnet is turned through a small angle Sip 
about a hne with direction cosines w, , 


8A,' = Sip 

by (49) of Chapter ii Hence, by (32), 
dW 

8w = sx: 


(34), 


= yilfilf' ^ {a. - 3 cos 6 


so that the couple about a hne parallel to the axis of is 




yMM' 


n |a. - 


„X, — X, 

3 cos 6 

r 


(35). 


Hydrostatics and Classical Hydrodynamics. The mtemal 
reaction in a fluid across an element of surface dS is a 
pressure pdS normal to that surface If the density is p, 
the bodily force per umt mass X, , and the velocity of the 
fluid at x^ IS w, , the acceleration of the fluid is found, by 
considering a small parallelepiped, to be given by 


dut 

P Ht 




( 1 ) 


If u, 18 given in the Eulerian way as a function of the 
coordinates x, and the time t, the operator djdt, giving the 
rate of change of any element associated with a given 
particle of the flmd, is equivalent to 


d _ a ^ 

dt~dt^ dx^ 


( 2 ) 


If we consider the circulation Q around any closed 
circuit C in the flmd and movmg with it, defined by 




I 


u,dx. 


( 3 ), 
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we have 


dQ _ 

dt 


f du, , 


L 


u. 


dt 


{dx^) 


The last integral is which always vanishes because 

when we move round the contour we come back to the 
same pomt, where the velocity has its original value Also 
if IS the gradient of a smgle-valued potential, as when 
the bodily forces are due to gravity (the commonest case), 

I Xidx^ IS the change of this potential round the contour and 


is zero Again, if is a function of p only, as in an incom- 
pressible liqmd or a gas at uniform temperature, the first 
mtegral vanishes and 



( 5 ) 


If then Q is ever zero around a circuit it remains so per- 
manently This 18 true if the flmd is mitially at rest and is 
set m motion by sohds movmg in it, and m various other 
cases of importance But the vamshmg of t 2 for all circuits 
IS the condition for the existence of a velocity potential <f) 
such that 


u, = 


d<i> 

dx. 


( 6 ). 


In this case we can rewrite the equations of motion in 
the form 


0W, 


du, I Sp „ 

dt 0Xfe p dx, ’ 

and multiplying by dx, and adding we have 


^^‘^dx +u —*dx - ---P 
dtdx.'^’^dxt ' pdx. 


dx, -t- X,dx, 


( 7 ) . 

( 8 ) , 


du, _ d^<f> _ duf, 

SXfc dx,dxj, dx. 


Binoe 


( 9 ). 



68 


CONTINUOUS SYSTEMS 


This shows that for all contemporaneous variations 

~ j ^ ^ 

where Xt = ^ (U); 

is the square of the resultant velocity q. The constant of 
integration is not necessarily the same at all instants and 
therefore may be a function of the time. Hence we have 
the Bemoulh integral 

+ + F{t) (12) 


The rate of change of mass withm a given small paral- 
lelepiped dx^dxj^dxg IS equal and opposite to the rate of 


outflow, hence 

^^dr = -^Jpu,)dr 


(13), 

and we have the equation of continuity 



II 

1 

(14). 

or 

dp 0M, 

dt ^ dx, 

(15). 


= A 


Vectors with given Divergence and Curl. We sometimes 
have to find a vector m, such that 

du, 

8x( 

dui dU( _ 
dx^ dxt, 

where A is a given scalar and tu,*. a given antisymmetncal 
tensor. We want particular integrals of these equations. 

d<j> 


( 1 ). 

( 2 ). 


Evidently if 
where cf> is any scalar, 




dxt 


(3), 


= A 


(4). 
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and a solution of (4) is 

where ^ is to be evaluated at a;, and A at ; 
dT = 

and r is the distance from to a:, 


where -F< is a vector such that 




we have 


' dx,dx^: ' 


Then also 

9% 9m, _ 9m* 

9i;~ aT*- '‘*“9x; 


by (9). 

Thus if 
we shall have 


_ ? M!? 

- 

a‘“j?’, 

d^F. 

— (a,j8„p— StpSms) g^ g^ 

- s _ s 

'"dXidx^ dx^dxj, 

9 9JF. 

9a:„ 9a:, 9a:, ^ 


-ft 111 


9m* 9m( _ 9m, _ 

9^;“ as; 9^; ~ 
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provided (9) is satisfied. But 

fff 

0a;„ 47rJJJ' 


“sill"- 35; (f)* 


since co„ is a function of alone and r a function of 
X, — Applying Green’s Theorem to all space except a 
small sphere about x, , we get 

dF„ 1 ^ ffLo.™ . 1 . rrr 1 , 


S-f’m 1 T 
2— = - 7- Lim 

cXm aTT 


,hi^dS+ Lim 
r Ait 


since Fn has the form of a gravitation potential and 
dF^jdx, that of a gravitational force But the first integral 
vamshes m the limit when the sphere becomes very small, 
and the integrand in the second is zero provided the com- 
ponents M, exist, for 

^ fix 

dx„ ?x„[dx, dxj '’■”‘dx,dx,„ 

If then we are given a scalar A and a vector such that 
its divergence is zero, a solution of (1) and (2) is 
dd dFn dF^ 

3 - + a“- 5 — (1^ > 

dx, dx^ dx^ 

where </> and F ^ are given by (6) and (12) and ikm are 
in cyclic order 

This analysis has two practical apphcations In hydro- 
dynamics M, IS the velocity and is twice the vorticity, 
denoted by 2^^ in Chapter ix Here the divergence of the 
velocity and the vorticity may be given through all space, 
and the velocity (17) satisfies the conditions If the actual 
velocity 18 v ^ , we may put 

V^ = U, + M,' (18), 


and then 


dx, ’ 


0M*' 5m/ 


- P- = 0 
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Thus u^' is the gradient of a scalar <f>' satisfying Laplace’s 
equation There is no such scalar that makes the velocity 
finite everywhere, including at an mfimte distance, except 
a hnear function of the coordmates Hence w/ is the same 
everywhere. 

If there are solid boundaries or free surfaces at a finite 
distance ii,' may not be constant If the region where 
vorticity 18 present does not extend to a boundary, there is 
no contribution to (14) from points outside this region, we 
therefore take ( 14) through a region large enough to contam 
the whole of the vorticity Then the surface integral in (15) 
must also be taken over the boundary of this region, but 
still vanishes, and (1") is still a solution But (17) may not 
satisfy the boundary conditions, and then we must add to 
w, an irrotational solution chosen to make the whole velocity 
satisfy them 

In electromagnetism m, may be the magnetic intensity 
and a)„, the electric current across unit surface in a plane 
of a;„ constant 

In many cases u, has no curl outside a limited region of 
very small cross-section This is often true m the motion of 
a real fluid, when the region may be called a vortex fila- 
ment, and m magnetism, when the region is a wire carrymg 
an electric current The former statement may be expressed 
also by saying that the motion is irrotational outside the 
vortex filament , the latter says that magnetic forces due to 
electnc currents have a potential In either case the mtegral 

= I u,dx, (20), 

taken around a closed circmt, is zero if the circmt can be 
filled up by a surface not cutting the filament or the wire, 
and has a constant value for aU circuits that cannot be so 
filled up The two conditions are mutually exclusive, and 
therefore the critical region itself must be a closed circuit 
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In each case A vanishes Then <^ = 0, and the component 
F„ is given by 



taken through the critical region. But if we consider the 
contribution from an element between two planes sepa- 
rated by df „ , and call the element of surface m a plane 
parallel to these dS, we have 

dr = di^dS (22), 

Ifco^ds = ju,di, = a (23), 


where the hne integral is taken around the boundary of 
the filament Hence 




(24), 


the integral bemg taken around the length of the filament 
Also 



(25), 


where ds is an element of length of the filament and a 
direction cosme 

In hydrodynamics is the circulation aroimd the fila- 
ment In electromagnetism the umt current is such that 
if it flows m a circle of radius 1 cm it produces magnetic 
force 277 at the centre. If we take the circle to be m the plane 
of Xg constant, with its centre at the origin, we have 


ds — dd , r = 1 , Xy= x^ — 0, 

= cos 6, ^2 == sm B, ^ 3 = 0 , = — sm 6, I 2 = cos 0, ls=0 
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and 


and 




ds = 0 for t = 


1 or 2 


= 2ir for 1=3 


Also for unit current Wa , the magnetic force, = 27r 
Hence, m this case, by (26), 

O = 4:TT, 

and m general, if the current is /, O = 4irl, and 



(26) 


(27) 


This may be transformed by Stokes’s Theorem mto an 
integral over a surface with the wire as its boundary , thus 


<■-11 


7 ^ ^ ^ri 

^tks ^3 


r r 

= - / 


i^iks 


(1) 

Vr/ 

jj Im — S.y S 


dS 


a* 


.( 1 ) 


dS 


= -ll'-Kv(-l)-5?k(l) 


dS 


11 


ai„ ai. (r) 


(28) 


But — 0 ^ '~ g ~ is the magnetic force at a:, due to 


a doublet of umt strength at f , with its axis in the direction 
Im Hence the force at any pomt is equivalent to that due 
to a distribution of doublets over the closmg surface, with 
intensity I per umt area and directed normally to the sur- 
face Such a distribution constitutes a magnetic shell of 
strength I. 
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It is customary to assert as a fundamental postulate that 
the magnetic force due to an electric current is equivalent 
to that of a magnetic shell, but I think this course un- 
desirable The magnetic shell does not exist in nature and 
direct experimental test is therefore impossible. Further, 
though it gives the same force, it does not give the same 
potential, the magnetic potential due to a magnetic shell 
IS a single-valued function with a finite discontmmty at the 
shell, while that due to an electric current is a cychc 
function with no discontmmty except at the wire itself 

Mutual Energy of Electric Circuit and Magnetic Field 
A magnetic pole m at a;, is under a force mu, , but 


mu. = — ml 


dx, 




and — is the magnetic force normal to the 


dS 


surface at i„, due to the pole Hence — j| 

IS the total flux across the closmg surface, or through the 
circuit, of the magnetic force due to the pole Also if F is 
the magnetic potential due to the current. 


u, = 


dv 

dx. 


(30), 


and therefore the correct forces are given by taking 


F = - / II NdS 


(31), 


where N is the normal magnetic force due to a unit mag- 
netic pole This can now be generahzed, we say that m 
general the mutual potential IF of a current and a set of 
magnetic poles is 


lF = -/||jVd>S 


(32), 
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where N is now the normal magnetic force due to the whole 
of the rest of the field This form may be further extended 
to express the mutual energy of two electric currents 
From this equation we can infer, as is done in the standard 


works of Jeans and Livens, that if N dS vanes with the 


time, the variation generates an mduced e m f 


dt 




NdS 


in the circmt. 

If we have two circuits carrymg currents I and J, their 
mutual influence is expressed by the statement that their 
mutual energy is of the form (32), where N is taken to be 
the force near due to the current J , pomts on the latter 
circuit can be taken to be given by x. Then 


W = IJM 


(33), 


where M = ^^l,u,dS (34), 


and w, IS the magnetic force at due to umt current in the 
second circuit But if F, is the vector potential due to such 
a current, 



(35), 


(36) 

= ^F„d^„ 

(37). 

taken around the circuit But 



F 


m 



(38), 


taken around the second circuit, and therefore 

dx^di„ 


M 


-II' 


(39), 


taken around both circuits This gives the required form 
for the coefficient of mutual induction of two circuits 


J C T 


5 
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^S7<- 

<3t > 

ISOTROPIC TENSORS 

A tensor is called isotropic if its components retain the 
same values however the axes are rotated We have 
already obtamed three examples, namely 6 , 4 , and 

There are no isotropic tensors of the first order For if 
u, was such a tensor, let us give the axes a small rotation 


expressed by the antisymmetrical tensor c,*. 
new system 

Then in the 

= (S.> - C,j) U,= U,- C,iUj 

and this can be equal to m, only if 

( 1 ), 

CjjUj = 0 

for all admissible values of the Cj, Thus 

^ 11^1 ^ 12 ^ ^13^3 “ 


( 2 ) 

^21^1 ^22^2 ^23 "^3 “ ^1 

^31^1 ^2^2 ^33^3 “ 


(3) 


But Cji = C 22 = Cg 3 = 0, while Cj^ , , Cjj are independent and 

equal and opposite to the components obtained by inter- 
changmg suffixes Hence (3) can be satisfied only if 


Uj — 112 — Wa — 0 ( 4 ), 

and therefore there is no isotropic tensor of the first order 
other than zero 

If u,t: IS an isotropic tensor of the second order, 

” ^tj) ~ ^kl) '“jl 

= '^ik ^kl^ti'^jl 

= Utk~ C^iUti ( 1 ) 

to the first order, for all values of ^ and k Hence 

CijUj^ + = 0 


( 2 ). 
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If i and k are unequal, take i = \,k = 2. Since Cii = Cgi— 0 
we have 

^ 12^22 " 1 “ ^ 13^82 “ 1 “ ^ 21^11 " I " ^ 23^13 ~ ^ (®)> 

and therefore 

^32 ~ ~ ® I ^U ” ^22 (^) 

By symmetry «,* is therefore Oiii ^ k, while = M 22 = “33 
If t and k are both 1, we have 

Ci 2'*^21 ^13^81 ^12^12 "f" ^13 ~ ^ (®)i 


which IS satisfied smee every term vamshes Hence the 
only isotropic tensor of order 2 is a scalar multiple of 8,;^ 
If u^^„ 13 an isotropic tensor of the third order, 

^ tkm ~ ^il) i^mn ^mn) '^jln (l)j 

and therefore, for all values of t, k, m, 

"I" ~ ^ (^) 

Take t = k = I Then 

®12^21m ^IS^Slm "i" •^12^12ni + 

“f" ~l~ ^m3^ll3 ~ 0 (3) 

Now put m — 2 so that 0^2 = 0 Then 

■^212 + Wj22 = Wijij 

'^312 + ^32 = '*hl2|’ (‘^) 

Wjl3 =0 ) 

Prom the last equation, and by symmetry, = 0 if two 
of 2 , k, m are equal and the third unequal Then by the 
first, M.tro IS also zero if all of t, I, m are equal, and the 
second shows that 

'^tkm ~ '^kz-m 

If in (3) we put m — every term vamshes, so that (3) 
holds 

Now in (2) if 2 , k, m are all different, is zero unless 
j = 2 , and then c„ = 0 Hence (2) holds It follows that the 


5-2 
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only isotropic tensors of order 3 are scalar multiples 

of 

If IS an isotropic tensor of order 4, we have, 

similarly, 

“1“ ^ ki‘^t jmn "t ^vs'^ikmi =0 (!)• 

There are only three possible values for *, m, p, and 
therefore at least two of them must be equal. We may con- 
sider separately the cases where (a) two are equal and the 
other two unequal, (6) three equal, (c) two equal and the 
other two equal, {d) all four equal. 

In case {a), take t = k=\,m=2, p=Z Then 

Oia'^aijs + Ci3 113128 + C 12 II 12 S 3 + 013111323 

-f- CiiUiiis -i- C23II1133 C3i11]i2i -4- C32ll]i22 = 0 (2). 

Hence, by the antisymmetncal property of c,* , 

%123 + %*23 ~ %U3 “ 

113123 + 11l323 — 'llll2l =0) 

lllisa — II1122 = 0 ( 4 ) 

Other instances of case (a) can be obtained by inter- 
changmg suffixes that are not already equal, and by tummg 
the axes so as to bung 3 mto the position of 1, 1 into that 
of 2, and 2 mto that of 3. Thus (4) gives 

■UllSS = ‘‘hm = II2233 = 112231 = “8322 = II3311 (5) 

And also 

lll313 = “1212 = «2323 = 112121 = 113232 = WsiSl (6), 
Hsus = ^2112 = 113223 = '^ 120 . 112332 = HuSl (7) 

In case (b), take 2 = k~m=l, p = 2 

C12II2112 -f Ci3 113112 + O12II1212 -(- O13II1312 

■b O12II1122 + O13II1132 -f- Cjglliui C 23 lliii 3 = 0 ( 8 ). 

The last term shows that 


iijuj — 0 


( 9 ), 
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and therefore, by interchange of suffixes, all components of 
class (6) are zero Also, from the coefficient of , 

'^112 '*^312 + '^hl32 — (f®)- 

But in (3) the last term vanishes and we infer 



^2123 ^223 — ^ 

(11). 

and therefore 


^312 "h W8U2 — 0 

(12), 

whence, by (10), 


'^132 — 0 

(13) 

Thus aU components 

of class ( 0 ) are also zero 



The coefficient of m (8) gives 


Wuil — ^2112 + '*^212 + *^U22 (14), 

SO that the components of class (d) are expressible in terms 
of the three types of class (c) 

No further information is got by transforming com- 
ponents of classes (c) and (d) Thus if i = fc=l,TO = p=2, 
replacing i or Ic by j will give a zero component unless j is 
equal to 1 , and then the factor c„ or c^, is zero, and the 
relation holds automatically Similar considerations apply 
if all of 2 , k, m, p are equal 

We may denote the components of type (5) by A, those 
of type (0) hy fx + V, and those of type (7) by g — v Then 
(14) gives 

■Wjjjj = 222222 ~ *^3333 ~ ^ + 2^ (15) 

There appear therefore to be three independent isotropic 
tensors of order 4, obtained by taking each of A, g, i/ m 
turn equal to 1 and the others to zero 

In the A tensor, — 1 if 2 ' = fc and m = p, and in all 
other cases is zero It is therefore equivalent to 
which is obviously a tensor of order 4, being the product of 
two tensors of order 2 

In the g tensor, = 1 if 2 = m, fc = p, or if 2 = p. 
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k — m, and i ^ k If also i = k, the component is 2 Other 
components are zero This can be written 

~ ^ip^km (I^)> 

and IS obviously a tensor of order 4 

In the V tensor, Ut^j, = liii = m, k = p, and = — 1 if 
I = p, k = m, and in all other cases is zero If also i = k, 
'^tkmp IS zero In this case, therefore, 

^ifcmB “ ^tm^kp ^tp^km 

This can also be written 

'^tkmp ~ ^tjk^mjp (^^)» 

for if t = 1, i = 3, = 0 unless = 2 and then = 1. But 

then e^jp =1 p = 3, — 1 ifm=3, ^)=1, and 

otherwise = 0 Thus 

■*^313 = ^SSl = — 1 (19)) 

With corresponding values for the other components. 
Evidently (17) and (18) represent a tensor of order 4 It 
has already appeared in Chapters I and vi 

The general isotropic tensor of order 4 is therefore 

+ fi + 8,b 8*,„) + v(S,„8tp- 8,3, 8t„) (20), 

where A, p, v are scalars. 


EXAMPLE 


Prove that 


“^tk ” '^mp '^pm* 

(^im^kp ~ '^mp ^pm 



CHAPTER VIII 


ELASTICITY 


In an elastic solid, as m a fluid, the distance between any 
two particles of the body usually vanes with the time. 
The body, however, has an eqmhbrium configuration that 
could persist if the external forces were zero or constant 
We may take this as a standard of reference. If a 
particle actually at x, would be at z, — u, m the standard 
configuration, we caU. w, the displacement at x, , in practice 
the squares of the u, can usually be neglected Evidently 
w, 18 a vector 

At a pomt X, + y, , where y, is small, the displacement is 
u, + V,, where 


du, 
= ax; 


ifiiK ^tk) Vk 


( 1 ) 

(2), 


where e,* and are the symmetrical and antisymmetncal 
tensors 


~ ^ vai, ax J ’ ^ V ax, ax J 


( 3 ) 


If 18 zero at x, , the displacement has the same form 
as that due to a general displacement u, together with a 
rotation expressed by the tensor Also, if e,^ is every- 
where zero, 

= 1 -- _ —A 

dx„ ^ dx„ \ax, dxj 

- 1 A. (A!^ 4 - i A / ^ -L 

^ ax, \ax„ axj^/ * dx^. \Sxn ax, / 

?^!!! — Q 

~ ~dxi axfc 


(4) 
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Hence the are constants and the rotation is the same 
everywhere. The vamshing of the e,^ is therefore the neces- 
sary and sufficient condition for a pure rotation. 

Now consider the part of v, due to the e,* . If we consider 
the quadnc surface 

6^ky^yk = (6), 

where ris a constant chosen so as to make the surface pass 
through y,, the normal at has direction cosines pro- 
portional to e,ky)c and therefore to The displacement 
due to e,t is therefore parallel to the normal at to this 
quadnc This comcides with the direction of the radius 
vector to y, if y, is on a prmcipal axis There are therefore 
three directions such that the relative displacement due 
to the IS in the direction of y, , and these directions are 
mutually perpendicular If we take new axes fj,, 
along them, the quadnc reduces to 

en'fi^ + + eaa'^s* = (6), 

and all terms e,i' with y # 1 are zero The displacement m 
the direction is now so that all distances in this 

direction are mcreased in the ratio 1 to 1 + e^'. The dis- 
placement due to the e,*. is then the resultant of three 
homogeneous strains parallel to three orthogonal axes 
We see therefore that the displacement m any small 
neighbourhood can be represented as the combmation of 
a rotation with three extensions at right angles The latter 
express the changes of size and shape of an element of the 
sohd. For this reason is called the strain tensor and 
the rotation tensor Evidently e,*. has six independent 
components For we can have 

= eyi, v^ = 0 (7), 

making all the e,*. zero except e^i = e, similarly 622 and 633 
may exist mdependently of the others Also, if 

ill = 0, = ez. Vs = ey 


( 8 ), 
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all the e,4 vanish except e^= e. Similarly e^i and Cja 
can be assigned independently. 

In an elastic solid the internal force across an element of 
surface is in general mclined to the surface If the area of 
the element is dS, the force across it must be specified by 
three components of the form PmdS parallel to the axes; 
where n may be regarded as indicating the normal to the 
surface. If we consider a small parallelepiped with edges 
dxi , dx2 , dxg , centred at {x^ , x^ , x^, the force across the face 
of area dx^dxg centred at + ^x^ is (pn, Pi2, Pis) dx2dx2, 
where the p^ are evaluated at x^ + \dxi . The force on the 
opposite face is — (p^, p^^, Pis) dx^dx^ evaluated at 

Xi — \dxi, and the total is dajjdxjdxa. In general the 

force m the direction due to the stress across the faces 

of iBfc constant is dx^dxidx^, and we take account of 

CXjt 

all faces by usmg the summation convention If the 
acceleration of the element has components /, and the 
density is p, the mass is pcteidxgdx,, while if the bodily 
force actmg has components X, per unit mass, the equa- 
tions of motion are 

The system of quantities p^^ constitutes a symmetrical 
tensor of the second order To prove this, consider first a 
plane whose normal has direction cosmos a „ , intersecting 
lines through x, parallel to the coordinate axes at short 
distances from x, , thus a small tetrahedron is formed, 
whose sides are of order I, say Let the area of the slopmg 
face be dS , then those of the others are a^jdS Consider 
now the rate of change of momentum of the matter within 
this tetrahedron. Evidently the contributions from the ac- 
celeration and X, are of the order of the volume, that is, 
of P. The force across dS has magmtude p„dS That across 
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the face of constant is — times the area of the 
face, that is, — p,ci^^id8. But dS is of order Hence we 
have 

{p,^-a^PM)Om = 0{l^) ( 10 ), 

and hence if Z is mdefimtely small we have for the stress 
across a plane normal to at x, , 

P,t = a^iPM ( 11 )- 

Now consider three jierpendicular di rections with direction 
cosmes {j = 1, 2, 3) The force per unit area across a 
plane perpendicular to one of these axes, in the direction 
of X, , IS given by p,, . Resolvmg this along the direction of 
one of the new axes Xj we have, therefore. 


Pii = = atjaiciPtic (12) 

by interchanging i and k, which is precisely the law of 
transformation of a second order tensor 

Consider again a small parallelepiped centred at x, , with 
edges parallel to the coordinate axes, and form the equa- 
tion of moments about its centre The contnbutions from 
fk and are of order I* at most, where the edges have 
lengths of order I The moment about an axis parallel to Xs 
of the stress in the face x^ = constant and parallel to x^ 
IS the product of p^ mto the area of the face and the dis- 
tance of the face from the centre, that is, to order Z®, 
^p^dx-^^dx^dx^. The opposite face makes an equal contribu- 
tion The stress pu in the face Xj = constant and parallel 
to Xi tends to turn in the opposite direction We have, 
therefore, 

(Pii - P12) dxidx^dxs = 0 {I*) (13), 


and therefore when we make Z indefinitely small we must 
have 


Pzi ~ Pvi 

and m general 

Pik ~ Pkt 

SO that the tensor p,/c is symmetrical. 


(14) , 

(15) . 
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Now consider the energy interchange between the small 
parallelepiped and its surroundmgs The stresses across the 
face of area dx^dx^ centred at x^ + ^dx^ are domg work on 
the element at a rate dx^dx^, and the contribution 

from the two opposite faces is 


dxi 


(Pti %) dx^ dx^ . 


Thus m all the stresses are domg work at a rate 
~{PuU.)dr 


The external forces are doing work at a rate pX,n,dr The 
kmetic energy of the element is ^pu^^dr, and is increasing 
at a rate pu^f^dr (We consider the actual specimen of 
matter occupymg the element of volume dr at time t , thus 
its mass 18 pdr, and if we keep to the same piece of matter 
at time t -\- dt the mass is unaltered. If we considered the 
variation of energy within a given element of volume we 
should have to allow for the variation of p and the fact that 
the matter moving out of the element is takmg its energy 
with it ) The rate of performance of work on the element 
therefore exceeds the rate of increase of kinetic energy by 



pXilli pu 



dr 




dr. 


by the equations of motion. 



dr 


(16) 


This work is stored up as internal energy of the element of 
volume Evidently from its form it is a scalar 

In any elastic sohd the internal energy is a defimte 
function of the state of the sohd In any change of state 
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the increase of internal energy therefore depends only on 
the initial and hnal states and not on the method of passage 
from one to the other Now we have seen that six of the 
e,i are mdependent, and if the element acquires displace- 
ments S-ii, in time the correspondmg morease of internal 

energy is 8 There is an apparent asymmetry 

according as i and k are equal or unequal Thus p^dr has 
coefficient S {du^ldx^) = Scu, but i = 1, A; = 2 contnbutes 
P218 (Swi/Saia), and i = 2, k = \ contnbutes PiaS (du^dx^), 
the two together giving 2pj^^8e^^. But this is the same as 
Pn 8^12 -f Pii Sfiji , and the whole contribution from the 
changes of stram is p^^Be^^dr Also during the process an 
amount of heat SQdr may be absorbed If then Edr is the 
internal energy of the element, 

BE — p,ic^,k + 8Q (17). 

Since E ve, & defimte function of the state of the system, 
and SIX of the Se,*. are independent and determine the 
other three, BE depends on the changes of the temperature 
and of the six mdependent e.^ and has a defimte value in 
whatever order these changes come about. But 


(dE\ __ 

(dE\ 

\0Cn/sQ=o 

\dei2'iQ=o 


2Pl2 


(18) 


If the absolute temperature is 6, and a certain amount of 
heat BQ is absorbed without change of any hnear dimen- 
sion, the nse of temperature is related to BQ by the rule 


BQ = pcBd (19), 


where c is the specific heat at constant strain If there is 
also a change of stram, smce BQ and 80 are scalars, we must 
have 

SQ = g.kSe,^-t pcSd (20), 

where the constitute a tensor of the second order. But BE 
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and 8Qj9 are perfect differentials Hence if we replace the 
SIX independent e.^ by e, we can wnte 




SE = {Pr + qr) SCr + pc86 

(21), 




(22). 

and 




+ ?.) = 

0 0 0 
= 0 ^^ (i>r + 3r). tPr + Qr) = 00 ^ 

(23), 


a 

3e, 

(9s\ 3 /gA a /q,\ a /pc\ 

\9) ~de,\9)’ d9\9) der\9 ) 

(24) 


It follows at once that if 0 is kept constant 2p,8e, and 
Sg'rSc, are perfect differentials Also 


a 

39 

{Pr + ?r) = 

9-^(^] 

der\9) 

= 0 - 
39 \t 

1 

II 

0 

(25), 

and therefore 


O^Pr 

^39 






9'r = - 



(26). 

If 80 = 

= 0 we can wnte 







Sc, 

,= Sir 



(27), 

where 







21F = 

Cq + 2c,er 

^rs^r^s 

+ 0(e*) 


(28) 


The Co, c,, Cr, may be functions of 9 Then if we retam 
only terms m IF up to order e*, 

Pr = c^ + Crse, (29). 

The Cf represent the stresses that would remain if the 
strams were removed without change of temperature 
In most practical cases the original state is one of uniform 
temperature and no stress, so that = 0 If there is a rise 
of temperature 9' under no stress, an element will acquire 
displacements 

v,= a,k9'yk (30), 
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where ai^ is a second order tensor expressing the thermal 
expansion Thus 

= i (“li + ®lbt) (31)> 

where /3,t is a symmetrical tensor, and 

Ptk ~ ^ik “h ^tk^mj ) (^^)) 

where „j, is a fourth order tensor 
But by hypothesis this deformation takes place under no 
stress and therefore p.f. ~ 0 This determmes 0 ,^, and our 
formula for the stress is 

Ptk = mj. (^mp — ^mp^') (33) 

The coefficient is the coefficient of in IF. 

Since there are six mdependent e, there are twenty-one 
possible terms m a quadratic form CrsCrCs, and therefore 
there are twenty -one coefficients c,* „„ They clearly form 
a tensor of order 4, such a tensor m general would have 
eighty-one components, but this satisfies the symmetry 
relations that it is unaltered if we interchange i and k, or 
m and p, or t and ^ together with m and p together 
From (26) and (33), 

0 0 

Qik ~ ^^mp ^0 i^ik^mp) 4 " ^ {^ik,mp^mp^ ) 

0 

~ ^^mp'^Q i^tk^mp) 4“ ^^ik^mp^mp (34) 

if 6' IS small The second term does not involve the e,*, 
the first IS small of the first order in the e.i. 

Many sohds are isotropic, that is, they have the same 
properties m all directions This applies to vitreous (glassy) 
sohds and to mixtures of crystals oriented at random In 
that case a uniform rise of temperature in an element 
gives an equal expansion m all directions and 

= aS’y, 


(36) 
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Simply ; then a is the coefficient of hnear expansion and 

^tk — ^^Ik (36) 

The second order terms m constitute a scalar, and 
we have 


Sp.k _ . 

o '^tk,mp 


(37), 


a tensor of order 4 If it is isotropic it must be of the 
form (20) of Chapter vn Then the hnear terms in give 

Ptk ~ ^ik mp^mp (38) 

= -f- fX {^tm^kp “b ^ip^km) ^mp 


"b ^ (3i)n ^kp ^km) ^mp 

~ "i" P "b ^zp^kp) ^ i^tm^mk ^zp^kp) 

= A B^k^mm p i^zk ■(“ ^kz) 4” *' i^zk ®li) 

= XBik^znm ~t" ^p^zk (39)> 


the last term vanishing since e,^ is symmetrical 
Then 


^W2 = Pzkezk (40) 

= Ae„emro + ^p^ik^zk 

= A A^+ 2,u (fill* + ^22^+ 633*+ 2e23^+ 2egi^+ 2ej2^) (41), 
where A = e„ = du^jdx^ (42) . 

The scalars A and p represent properties of the material. 
Both are positive 
We can also write 


21^2 = (A + 2p) A® + 4/Li (e23® + 

®22®33 ®33®11 “ ®11®22) (43). 

This appears to differ from the form in Love’s Elasticity, 
1906, p 100, but the present e,*, differ from Love’s stram 
components My en is the same as his, namely du^jdXk , but 
my fijs IS only half his, so that his assemblage of strain 
components is not a tensor 
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If all the e„ were equal to one another and therefore to 
^A, we should have 

2W,^{X+ Im) (44). 

= (^ + f/i) A (» = &), p,k=0 {i¥^lc) (45). 

In general we write 

A + §/. = A: (46), 

and caU k the bulk-modulus. 


2 W^= (A 

-f- 

f^)A^ 








-f 

2/1 (^11^ 

+ 

^22^ 

-b 

633* -F 

2p 2 

■^^23 

+ 2 e,,^ + 2 e,/) 


- 

1/1 (^11^ 

- 1 - 

^22^ 

- 1 - 

Caa^-F 

2^22 ®83 ^^ 33^11 ”S~ ^^11^2) 

= (A 

-f 

<1 








+ 

t/l 

+ 

fi 2 
^22 

-F 

£33*- 

622^33 

— 653611 — 611622) 


-F 

4/1 (Caa® 

•f 

fi 2 
^31 

"F 




= (A 

-f 

f/i)A^ 








- 1 - 

f/l {(^22 

- 


' + 

(^ 

eix)^ 

+ (®ii “ ^2)^ 


-t- -I- eegi* + 6ei2*} (47). 

The coefficient of n vamshes if and only if the strain is 
a symmetrical expansion, and may therefore be called the 
distortional strain energy 

If we allow for variations in temperature, 

Pik = A (A - 3ad') + 2fi (e,* - a^'S,*) 

= {AA - (3A + 2p) ad'} + 2^e., (48) 

and g-a = - ^ (49) 

Every term in is of the first order in the displacements , 
but dd'/dd — 1 and therefore gives rise to a constant term. 
This term is 


9ik ~ ^ ■3ka8fic 


(50) 
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If SQ = 0 , 80 that no heat is lost or gamed by conduction, 
pcBd= - g'„;8e,fc= - 

= - SkadSA (61), 

and therefore, if the stram takes place adiabatically, 

^kad . 


and p,, = (^A + 8., A + 2pe,, (63) 

Thus in an adiabatic disturbance the constant A is in- 
creased above its value for a standard disturbance to A', 


where 


A' = A-F 


9k^a^e 


while fj. 18 unaltered The bulk-modulus k is therefore also 
mcreased to k', where 

r = + (55). 

pc 

In a simple thermal expansion at zero stress the absorp- 
tion of heat 8Q is equal to pCj,S0, where Cp is called the 
specific heat at zero stress, and is the specific heat mea- 
sured m ordmary experiments Then 

pCp86 = 8Q = pc 80 + g'.i8e,fc 

= pcS0 + 3ka0Sji^8ei^: 

= pc80 -|- 3k a0 3a80. 
rr,, 9ka^0\ kc 


Cp = c 1 -f 


kc 

rk' 


The equations of motion at constant temperature, if the 
properties A and /x are umform, can be written 

p!% = + 2/xe.J + pX, (57) 

9 . 9 /duu . 




~ ST + m) ^} + + P^t 


J C T 


6 
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If the changes are adiabatic A must be replaced by 
A'. 

If there is any heat conduction, the absorption of heat 
per umt tune by the element of volume dr is 



where K is the thermal conductivity Then the equation 
of heat conduction is 


(59) 

If we write 

e,j. = ciS,id' + e,; ' 

(60), 

so that e,j' is the stram due to the stresses, 


9ka^0~+3kae^' 

(61), 

and the equation becomes 


+ 3kae~==KV^e 

(62) 

The element dr originally had volume 


0 (Xj llj , Xj U 2 , Xj M 3 ) 

0 ) ^2 > ^ 3 ) 


"{‘-(S-zaM”:) 

dr (63), 

SO that its density was p (1 + A) Hence 


dA 1 dp 

dt p dt 

(64) 


This is the equation of continuity 
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HYDRODYNAMICS 

In comparison with a typical elastic sohd, a real fluid shows 
a great resemblance and a fundamental difference The 
force per umt area across an element of surface parallel to 
a coordinate plane again constitutes a symmetrical tensor 
of order 2, and for the same reasons The equations of 
motion are still (9) of Chapter vin, and also the rate of 
performance of work on an element of volume dr still has 
the form (16). The difference is that the mternal energy in 
a flmd does not depend directly on how much it has been 
deformed However the fluid is moved about and stirred 
up, provided it returns to its mitial position, density, and 
temperature, the imtial and final internal energies are 
equal The deformation, however great it may be, makes 
no oontnbution , the stresses do work on each element, and 
thereby supply energy, but this is removed in restoring the 
onginal temperature If energy of deformation existed in 
a fluid, all particles of it would have a tendency to return 
spontaneously to their onginal positions when stresses are 
removed, and they have none Accordingly, while the rate 
of performance of work on the element of volume dr m time 
dl IS still picjdujdxicdr, where the symbols have the same 
meamngs as in elasticity, we can no longer assert from this the 
equation (18) of Chapter vin, because the change of internal 
energy is not determmate when the changes of the are 
given The e,^ may be as great as we hke, but the energy 
does not increase indefimtely apart from changes m density 
and temperature ; and the flmd moves in the same way under 
the same external forces whatever its previous history. We 
may say that an elastic sohd has a memory , a flmd has none 

6-2 
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The stresses in a fluid are related, not to the total de- 
formation, but to the rate of increase of deformation. In 
the former notation these have components 


Sc.j 

"dt 


or 


1 j. 


The velocities now appear, mstead of the displacements 
from an imtial configuration, and we now denote the 
velocities, instead of the displacements, by u, . 

We also wnte 


so that e.fc now denotes the raie of increase of strain and 
the local angular velocity, usually called the vorticity 
Now we say that p,^, is hnearly related to the new e,*., and 
therefore 


Ptk "b ^ik,mv^mp (^)> 

where P,* is a symmetrical tensor of order 2 that possibly 
does not vamsh with the e,* , and c.j, is a tensor of order 4 

Further, in a fluid at rest the stress is an isotropic tensor. 
Hence P,* is isotropic and must be of the form — pS,^, 
where p is a scalar, and c,^_„j, must be of the form (20) 
of Chapter vn On carrying out the summation with regard 
to m and p the term in v disappears and we have 

Pik ~ P^tk "b (3)) 

where A and p are scalars, at present undetermined 
Evidently p and A are not both necessary to preserve the 
form, and we may introduce the further convention 


P,t = - 3p (4), 

so that — p IS the mean of the three p,*. with equal suffixes 
This gives 

- 3p = - 3p -f- 3Ae„„ + 2pe„ (5), 

and therefore A = — §p. (6). 

If we consider the internal energy as a function of the 
density p and temperature d alone, we may consider the 
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energy change in a symmetrical expansion. If the density 
mcreases by Sp, there is a contraction m all dimensions in 
the ratio ^8p/p and the stresses do work 

At the same time there may an absorption or generation 
of heat , then the energy change is 

+ (8), 

P 

and the heat absorbed may be written 

8Q = Mhp + pc89 (9), 

where M is unknown and c is the specific heat at constant 

volume. Then \ 

+ JUj 8p + pc8d (10), 

( 11 ), 

and the condition that these quantities are perfect differen- 
tials gives 

It appears also that p must be a function of p and 0 
alone, for a given material , it does not involve dp/dt * We 
may call p the pressure 


* This amounts to saying that there is no dissipation of energy in a 
symmetrical compression or expansion This is true in a gas on the older 
kinetic theory , but Enskog has shown (Kungl Svenska Akad Handhngar, 
63, no 4, 1922, p 18) that p can with greater accuracy be given by 


where y is a “second ” coefficient of viscosity But 17 //X is only of the order 
of the square of the ratio of the volume of the molecules themselves to the 
whole volume of the gas 

In a liquid the coefficient of dp/dt, if it exists, is within the experimental 
error, an analogous statement is true for imperfectly elastic solids 
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The stress components may now be written 

Pik = - (p + 8.* + 2/xe,fc (14), 

where p is a function of the density and temperature and 
fj. expresses a property of the flmd. We call n the coefficient 
of viscosity For a uniform flmd we have the equations of 
motion 

^ {“ ^)} 

9 9 

= - 0^ (i’ + 0— e** 

= (15), 

where we now write the scalar 


= A 


( 16 ) 


In time dt the outflow from a volume element dkc^ dx^dx^ 
is {0 (pw,)/0a;,} drdt The mass within the element therefore 
decreases at a rate {0 (pUt)ldx,} dr, and we have the equation 
of continuity 

dp 0 


or 


d 


dt dx, 
dp du, . 

jr-Pdxr-^^ 


( 17 ), 


( 18 ) 


Here r denotes differentiation with regard to the time 
at 


dx 

following a given particle of the flmd, so that 
given by 


are 


dxt 

dt 


= u. 


and 


d 

dt 


0 

dt 


(19) , 

( 20 ) . 


dx, 8 _ 0 0 

dt dx, dt dx. 

In most works on hydrodynamics this operator is denoted 
by D/Dt, but I see no reason for departmg from the usual 
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notation for total differentials, since this particular one 
18 the only total differential that occurs in hydrodynamics. 
The acceleration components are given similarly by 


/.= 


dUt 

dt 


du. 


du. 


dt ~dx^ 


( 21 ). 


Consider now the circulation O about a closed contour. 
We have 


dO. 

dt 


The integral I always vanishes, | vamshes 

for bodily forces denved from a potential Then 


■ ~ L P al L p a*; 

ipA) + j J A 2f.. + 2 g-*)} .fe. 

- - 1 p ■* <? + 5'^) -IJi, 


+ I - 
3c P 


2^ 

dXi 


d'Uit+ - p.d^ ( 23 ) 
3c P 


If /I = 0 this reduces to the circulation theorem of classical 
hydrodynamics In many problems of real flmds ju, is 
small and constant, and A small Then the mtegrals 

[ -rf(|yu.A), f - ^dui, and [ -/iirfA are zero or products 
J Cp J C P J C P 

of two small quantities, and can be ignored. - is zero 

3 cP 

if p 18 a function of p, which is true in many cases, though 
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exceptions arise when the temperature or the composition 
vanes from place to place. In the commonest case, to 
considerable accuracy, 

and circulation anses only from vanations of vorticity in 
the neighbourhood of the contour It follows that in a 
flmd originally at rest or in irrotational motion circulation 
can anse only through the diffusion of vorticity inwards 
from a boundary 

The work done on the element dr per unit time exceeds 


the rate of increase of kinetic energy by 



Pu = Pti^e.^dr 

(25), 

and 

Ptk^tk = — (P + IpA) + ^P^tk^ik 



= - (P + |p^) ^ 

(26) 

But 


(27), 


and if p 18 a function of p (including a constant as a par- 
ticular case) this is a differential with regard to the time 
and yields on the whole zero if the original density is ever 
recovered This term may therefore be considered to give 
the increase of internal energy due to compression The 
remainder may be written 

O = - §pA2 + 2pe,je,s 

~ fp {(® 2 Z ~ ^33)^ + ®ll)^ + (®11 ^22)*} 

+ (6232 + 631^ + 612*) ( 28 ) 

Thus d> IS analogous in form to the distortional strain 
energy of elasticity ((56) of Chapter vra) It is essentially 
positive, and therefore represents work done on the flmd and 
continually stored as internal energy. In (8), therefore, if the 
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change is an actual one, the term (27) contributes thepSp/p, 
while (28) contributes to 8Q So long as the imtial and final 
states are given, it is immaterial whether 8Q represents 
heat conducted into the element, absorbed from radiation, 
or generated chemically within it, or mechamcal energy 
dissipated into heat by viscosity. 

We notice that <& can vamsh if, and only if, 

^11 ~ ^2 ~ 

623 = 631 = C12 ~ (2^)> 

so that the deformation represents a symmetrical expres- 
sion or contraction. 

If we consider any fimte volume, the rate of dissipation 
within it IS 


III ^dr = III (- fyxA* + dr (30) 

But = + (31), 

and 



(32), 


where is a direction cosine of the outward normal to the 
boundary But 



=^jj2l,uJ,,dS+ijj^-^'dS 

= II 2e,,^uMmdS + i II dS (33), 
where d/dn denotes differentiation along the outward 
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normal, q is the resultant velocity, and is the vorticity 
vector associated -with Also 



(34), 

and m all 

jjj Odr = dfS + 4/i |j 


+ 2^ III (- JA* + ^ dr (35). 


This form is useful m such a problem as that of waves 
on deep water, where the viscosity is small, if it were 
absent we should have a permanent oscillation in a normal 
mode The vorticity is neghgible except within a distance 
from the bottom of order (v/y)K where = vp and 2iTjy is 
the period of the motion Then in the second integral 
IS zero at the free surface and «, is zero at the bottom, so 
that this integral vamshes A is neghgible everywhere 
Near the bottom contains a factor proportional to the 
velocity that would exist there on the classical theory, and 
thus can be made mdefimtely small for deep water. Hence 
the important term in (35) is 

and even this vamshes at the bottom, so that it need only 
be estimated at the free surfaxie 

If a portion of the fluid is compressed without change of 
temperature, we define the bulk-modulus k by 


Idp _l 
p dp k 



and if it expands under change of temperature without 
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change of pressure we define a coefficient of volume 
expansion a by 

1 %--“ (*’> 

Then for small changes of temperature and pressure 




Also m a free expansion under constant pressure 
SQ = M8p + pcSd 




where the partial differentiation is to be earned out at 
constant density But this makes 



- ak 
dd - 

(40), 

and therefore 

SQ = p |c + a* ^ 89 

(41), 


so that Cj, , the specific heat at constant pressure, is given by 

. a^kd 


This form is analogous to that of (56) of Chapter vin, the 
present a being the coefficient of volume expansion, and 
the previous one that of linear expansion 
In an adiabatic change SQ = 0, then 

— akd^ + pc86= 0 (43), 

P 

and therefore ad8p = (pc + a^kd) 89 

/, a^k6\ akS „ 


Sp = ^: ^1 + 


-- bp, 
P 


pc I p 

a^k9\ 8p 
pc J p 


so that 


(44) 
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Thus the bulk-modulus for adiabatic changes is 
,, , . a^kd\ kc„ 


\ pc J c 


The equation of heat conduction needs to be modified 
to allow for heat generated internally In time St, per 
volume dr, we have 

= pcSd -|- akdASt (46), 

whence + nMA - O + A (iT |) (47) 


EXAMPLES 

1 Obtain the equations of motion in terms of the stress 
components by considering the momentum of a finite volume of 
any form and applying Green’s theorem 

2 Similarly obtain equation (25) 
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Acceleration, 24 
Angular momentum, 35 
Area, plane, as vector, 20 
Axes, change of, 1 
dynamical, 24 
rotatmg, 40 
Axis, central, 45 

Centroid, 33 
Change of axes, 1 
Circulation, 56, 87 
Couple, 43 
Curl, 14, 58 

D’Alembert’s prmciple, 32 
Dielectrics, 52 
Direction cosines, 1, 16 
Dissipation, 88 
Divergence, 8, 58 

Elasticity, 71 
Electromagnetism, 61 
Electrostatics, 60 
Energy, kmetic, 25 
potential, 26 

Equation of contmuity, 58, 82 
Force, 26 

Ime of action of, 42 

Gauss’s theorem, 49 
Gravitation, 49 
Green’s theorem, 48 

Hamilton’s principle, 28 
Heat conduction, 82 
Hydrodynamics, 83 
Hydrostatics and classical hydro- 
dynamics, 66 

Inertia, moments and products of, 37 
Isotropic solids, 78 

Kmetic energy, 25 


Lagrange's equations, 28 

Magnetism, 53 
Moment of force, 43 
of vector, 36 

Null planes and pomts, 46 

Parallel, 17 
Potential, 49 

Rigid bodies, 29 
Rotatmg axes, 40 
Rotation, 21, 30 

Scalar product, 7 
Scalars, 3 

Specific heat, 81, 91 
Stokes’s theorem, 48 
Stram tensor, 72 
Stress tensor, 73 
Suffixes, dummy, 5 
Summation convention, 3 

Tensors, 2 

antisymmetrical, 5, 8 
isotropic, 66 
symmetrical, 5 
Kky 6 

^jkm> 

Vector product, 7, 14 
Vectors, 2 

with given divergence and curl, 
58 

Velocity, 24 
Virtual work, 34 
Viscosity, 86 
Vorticity, 60 

Work function, 26, 50 
Work, virtual, 34 
Wrench, 45 
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